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1.0 Introduction 


METCAN (Metal Matrix Composite Analyzer) is a computer program developed at 
NASA Lewis Research Center (References 1-3) to simulate the high temperature 
nonlinear behavior of continuous fiber reinforced metal matrix composites. 

METCAN incorporates constituent material models along with composite 
micromechanical and macromechanical models to allow a comprehensive point 
analysis of the composite thermal and mechanical behavior. 

The following sections provide a step by step outline of the procedure necessary 
to run METCAN. Section 1.0 will begin with an overview of METCAN, followed by 
a sample METCAN problem, the various geometric conventions employed, the 
structure of the input/output files, a description of the micromechanical unit 
cell model and the nonlinear constitutive relationship used in METCAN, and a 
list of notation. Section 2.0 will demonstrate the preparation of the primary 
input file, Section 3.0 will show how to create the constituent databank, and 
Section 4.0 will explain the contents of the postprocessing files. 

The preparation of the following METCAN companion manuals is also planned 
for the near future: 

i) METCAN Demonstration Manual: This manual will contain various 
sample problems and details of the corresponding input and output files. 

ii) METCAN Theoretical Manual: This manual will detail the theories 
used in the development of METCAN. 

iii) METCAN Programmer's Manual: This manual will provide general 
programming information as well as describing the contents of the various 
subroutines that compose METCAN. 
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The potential user of METCAN is reminded that the program and the 
theoretical methodologies implemented in the code are in an evolutionary 
state* The present release of Version 4.0 is on a preliminary basis for 
evaluation purposes. Additional information regarding METCAN and the 
efforts to validate and verify the code can be found in References 4-8. 


1-3 



1.1 METCAN Overview 


High temperature metal matrix composites offer great potential for use in 
advanced aerospace structural applications. The realization of this goal 
however, requires concurrent developments in (1) a technology base for 
fabricating high temperature metal matrix composite structural components, (2) 
experimental techniques for measuring thermal and mechanical characteristics, 
and (3) computational methods to predict their behavior. In the development of 
high temperature metal matrix composites, it proves beneficial to initially 
simulate their behavior through computational methods. Apart from providing 
an initial assessment of the metal matrix composite, this method helps to 
minimize the costly and time consuming experimental effort that would otherwise 
be required. 

Recent research into computational methods for simulating the nonlinear behavior 
of high temperature metal matrix composites at NASA Lewis Research Center has 
led to the development of the METCAN (Metal Matrix Composite Analyzer) computer 
code. METCAN treats material nonlinearity at the constituent (fiber, matrix, 
and interphase) level, where the behavior of each constituent is modelled using 
a time-temperature-stress dependence. The composite properties are synthesized 
from the constituent instantaneous properties by making use of composite 
micromechanics and composite macromechanics models. Factors which affect the 
behavior of the composite properties include the fabrication process variables, 
the in-situ fiber and matrix properties, the bonding between the fiber and 
matrix, and/or the properties of the interphase between the fiber and matrix. 

The METCAN simulation is performed as a point-wise analysis and produces 
composite properties which can be incorporated into a finite element code to 
perform a global structural analysis. After the global structural analysis is 
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performed, METCAN decomposes the composite properties back into the localized 
response at the various levels of the simulation. At this point the 
constituent properties are updated and the next iteration in the analysis is 
initiated. This cyclic procedure is referred to as the integrated approach to 
metal matrix composite analysis and is depicted in Figure 1.1-1. 
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Figure 1 .1-1— Integrated approach to metal-matrix composite analysis. 
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Figure 1.1-2 shows the modular structure of METCAN. In the development of 
METCAN, emphasis has been placed on maintaining a modular software structure and 
in providing a user friendly interface. The code features (1) a dynamic storage 
allocation scheme for efficient use of computer resources, (2) a resident 
databank of constituent material properties, (3) user selected control of the 
printed output, (4) generation of postprocessing files from which selected 
output can be extracted and represented graphically, (5) an input file structure 
which provides a straightforward user interface, and (6) separate modules 
containing the failure criteria, the material model, the composite 
micromechanics analysis, and the laminate analysis which are incorporated into 
METCAN. 



Figure 1 .1 -2 — Modular structure METCAN. 
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1.2 An Initial METCAN Sample Case 


A sample case with a step by step description of the procedure required to 
prepare a METCAN input data file is provided below. This sample case is 
presented at this point to give the user an overview of a typical application of 
METCAN to model a problem before delving into specific details, as well as to 
provide a starting point from which the user can easily utilize METCAN. The 
problem to be modelled involves a two ply unidirectional metal matrix composite 
composed of SiC fibers and a Ti-15-3 matrix. The fiber volume ratio is 30% ,no 
voids are present, each ply has a thickness of 0.005 inches, and a perfect bond 
(i.e. no interface) is assumed between the fiber and matrix. The composite is 
processed from a temperature of 1200 F to room temperature and then subjected to 
a monotonic load in the longitudinal direction at room temperature, as shown in 
Figure 1.2-1. Results taken from the output containing the composite moduli at 
the end of processing and at the end of loading are also included to demonstrate 
one of the typical capabilities of METCAN. 
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Figure 1 .2-1 — Sample case loading history. 
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The individual steps along with the corresponding line numbers in the input data 
file required to model this problem are presented in Table 1.2-1, while the 
input data file formed from these steps is shown in Figure 1.2-2. The line 
numbers shown in the table and figure are included for ease of reference between 
the two and are not part of the input data file. Typically, any problem can be 
modelled in ten steps, where each step requires one or more lines of data. 


Line No. 

4- 






+ 



1 

SAMPLE CASE 

FOR A 

TWO PLY 

SiC/Ti-15-3 METAL 

MATRIX COMPOSITE 




2 

POST 

F 








3 

LDDIST 

F 








4 

PLY 

1 

1 

0.0 0.005 






5 

PLY 

2 

1 

0.0 0.005 






6 

MATCRD 

1 

.30 

0. S I CAT 115 






7 

CYCLES 

1.0 

1.0 







8 

PRINT 

ALL 








9 

PRINTOPT 

ALL 








10 

INTRFACE 

0 








11 

TMLOAD 

-2 








12 

0.0 

1000. 

35 







13 

1200. 

1200. 








14 

70. 

70. 








15 

0 . 

0 . 

0 . 

0 . 0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

16 

0 . 

0 . 

0 . 

0 . 0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

17 

1000. 

2000. 

25 







18 

70. 

70. 








19 

70. 

70. 








20 

0 . 

0 . 

0 . 

0 . 0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

21 

500. 

0 . 

0 . 

0 . 0 . 

0 . 

0 . 

0 . 

0 . 

0 . 


+ + — - — + + + + — + + — + 

Figure 1 .2-2 — Sample case Input data file. 
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Table 1.2-1: Procedure in Preparing an Input Data File 


Step No. Line No. 


Description of Data 


1 

2 


4 


5 


1 Choose a problem title. 

2 On the ’POST* record, select ’ F’ to indicate that the 
postprocessing files (used for graphical representation of 
the output) are not to be generated. 

3 On the 1 LDDIST* record, specify ’ F* to show that the load 
redistribution algorithm (a post failure model used to 
simulate the progression of damage in the composite) is not 
to be used. 

4 Use the ’PLY 1 record for the first ply. Specify the ply 
number as 1, the * MATCRD* record identification number as 1, 
the ply orientation angle as 0°, and the ply thickness as 
0.005 inches. 

5 Use the ’PLY* record again for the second ply. Choose the 
ply number as 2, the ' MATCRD* identification number as 1 
(indicating that both plies make use of the same * MATCRD’ 
record), the ply orientation angle as 0°, and the ply 
thickness as 0.005 inches. 

6 Specify on the MATCRD’ record the material identification 
number from the * PLY' records at 1, the fiber volume ratio 
as 0.30, the void volume ratio as 0.0, and the fiber and 
matrix code names used to retrieve the constituent material 
properties from the databank as SICA for the SiC fiber and 
TI15 for the Ti-15-3 matrix. 
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6 


7 


Use the ’CYCLES’ record to select the number of thermal and 


7 8 


8 9 


9 10 


10 11 


12 


13 


14 


15 


16 


17 


mechanical cycles as 1.0, indicating that no degradation in 
material properties due to cyclic loading is desired. 

Use the ’PRINT’ record and specify ’ALL’ to print all the 
possible output requests. 

On the ’PRINTOPT’ record, choose ’ALL’ to indicate that the 
output requests from the previous card are to be printed 
after each step in the analysis. 

Use the ’ INTRFACE’ record with a ’0’ to specify that an 
interface is not modelled. 

Use the ’TMLOAD’ record to designate the loading history. 
Selecting a ’-2’ indicates that a linear loading history 
split into two segments is to be used. The first segment 
simulates the processing of the composite, while the second 
segment models the application of the mechanical load. 

Set the time range for the first segment from 0 to 1000 
seconds and assign this segment 35 increments. 

Specify the temperature of plies 1 and 2 at the start of the 
first segment (0 seconds) as 1200 F. 

Specify the temperature of plies 1 and 2 at the end of the 
first segment (1000 seconds) as 70 F. 

Indicate that there are no mechanical loads at the start of 
the first segment by choosing ’0.0’ for all entries. 

Indicate that there are no mechanical loads at the end of 
the first segment by choosing ’0.0’ for all entries. 

Set the time range for the second segment from 1000 to 2000 
seconds and assign this segment 25 increments. 
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18 Specify the temperature of plies 1 and 2 at the start of the 
second segment (1000 seconds) as 70 F. 

19 Specify the temperature of plies 1 and 2 at the end of the 
second segment (2000 seconds) as 70 F. 

20 Indicate that there are no mechanical loads at the start of 
the second segment by choosing ’0.0* for all entries. 

21 Indicate that at the end of second segment there is a 
longitudinal force resultant of 500 lb/in applied. This 
means that the load increases monotonically from 0 at 1000 

seconds to 500 lb/in at 2000 seconds. 
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Once all ten steps have been completed and the input data file generated, the 
METCAN program can be executed. Results of the METCAN simulation will in turn 
be written to an output file. Typical results from the output file for the 
composite moduli at two points in the loading history are shown in Table 1.2-2. 
A more detailed version of this input data file along with two other input data 
files can be found in Section 2.13. 


Table 1.2-2 Typical Results for Composite Moduli Taken from METCAN Output 


Moduli (Msi) 

End of Processing 

End of Loading 

E 

Cll 

26.63 

24.97 

E 

C22 

16.66 

16.52 

E 

16.66 

16.52 

C33 



E 

C12 

6.61 

6.55 

E 

C23 

6.26 

6.23 

E 

C31 

6.61 

6.55 
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1.3 METCAN Geometric Conventions 


The global and local coordinate systems employed in METCAN and the 
conventions for ply ordering, measuring orientation angle, and loading are 
presented in the following sections. 
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1 . 3 . 1 


Global (Laminate) Coordinate System 


Figure 1 - 3 . 1 1 shows the global or laminate coordinate system for an arbitrary 
laminate. The convention is that the x-axis is orientated parallel to the 
direction of the fiber for a 0° ply, the y-axis is orientated perpendicular to 
this fiber direction, and the z-axis is orientated through the thickness 
dimension of the composite. The x-axis direction is also referred to as the 
longitudinal direction, while the y-axis and z-axis directions are called the 
in-plane transverse and out-of-plane transverse directions, respectively. 



Figure 1 .3.1 -1 METCAN global (laminate) coordinate system. 
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1.3.2 


Local (Material) Coordinate System 


The local or material coordinate system is depicted in Figure 1.3. 2-1. The 
1-axis is orientated parallel to the direction of the fiber, the 2-axis is 
orientated perpendicular to this fiber direction, and the 3-axis is orientated 
through the thickness dimension of the ply. 



Figure 1 .3.2-1 — M ETC AN local (material) coordinate system. 


1.3.3 


Ply Ordering Convention 


The ply ordering convention employed in METCAN is shown in Figure 1.3. 3-1. The 
first ply is placed at the plane z=0. Subsequent plies are stacked on top of 
the first ply in the positive z direction. 



Figure 1. 3.3-1— METCAN ply ordering convention. 
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1.3.4 


Orientation Angle Convention 


The orientation angle of a ply is taken as the angle between the local 
coordinate system and the global coordinate system, as shown in Figure 1.3. 4-1. 
The positive convention for angle 0 is measured counter-clockwise. 



X, Y : global coordinate axes 

1,2: local coordinate axes 

0 : angle between global and local axes 

Figure 1 .3.4-1 — M ETC AN ply orientation angle 
convention. 
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1.3.5 


Loading Conventions 


The positive sign conventions for the various loadings in METCAN are presented 

in Figure 1. 3.5-1. It should be noted that the METCAN input file assumes 

stresses and moments to be integrated resultant quantities in the usual sense 

from plate and shell theory. Dimensions for the force resultants are force per 

unit width, while moment resultants are moment per unit width. The in-plane 

force resultants are represented by N , N , and N , while out-of-plane force 

x y xy 

resultants are represented by Q and Q . Moment resultants are symbolized by 

xz yz j j 

M , and ^ , and the upper and lower surface pressures are denoted by P and 

x y xy u 

P^, respectively. 



Figure 1.3. 5-1 — METCAN positive loading conventions. 
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1.4 METCAN Input/Output File Structure 


The input/output file structure in METCAN is comprised of four different files 
as shown in Figure 1.4-1. Three out of the four files (the two input files and 
the primary output file) will be utilized for all analyses. The fourth file is 
actually a set of ten postprocessing files which can be optionally generated as 
defined by the user. The primary input file contains information regarding the 
ply material system, the loading history, and identifies the constituent 
properties to be retrieved from the second input file. The second input file is 
called the constituent databank and contains the room temperature constituent 
properties for various materials. Data from the two input files is used to 
initiate the METCAN analysis. Results from the analysis are in turn generated 
as the primary output file. The primary output file can be tailored by the user 
(as defined in the primary input file) to produce various combinations of the 
available output requests. The available output requests include constituent 
mechanical and thermal properties, micromechanical stresses and strains in the 
three subregions of the unit cell for each ply, ply mechanical and thermal 
properties, and composite properties. Another group of ten output files, the 
postprocessing files, can also be optionally generated (again defined through 
the primary input file). The postprocessing files contain the same information 
found in the primary output file, but in a format more amenable to manipulation 
for graphical representation. The following sections give a brief description 
of each input and output file. 
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Figure 1 .4-1 — M ETC NAN input/output file structure. 
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1.4.1 


Primary Input File 


The primary input file contains data used to initiate the METCAN analysis 
The file performs the following four functions: 


1) Specify the composite geometry. Variables such as ply thickness, 

orientation angle, fiber volume ratio, and interphase thickness 
are prescribed. 


2) Specify the loading history. The fabrication conditions. 

thermal and mechanical loadings, and thermal and mechanical cyclic 
loadings can be detailed. 


3) 


Identify the constituent materials. The code names of the fiber, 
matrix, and interphase to be retrieved from the constituent 
databank are specified. 


4) 


Select output information generation. Allows the user to 
the primary output file to the user's requirements and to 
the generation of the ten postprocessing files. 


tailor 

control 


These functions of the primary input file are depicted 
and the specific data records necessary to prepare the 
will be addressed in Section 2.0. 


in Figure and 
primary input 


1.4. 1-1 
file 
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Figure 1. 4.1-1 — Functions of the primary input file. 


1-23 














1.4.2 


Constituent Databank File 


METCAN utilizes a constituent databank containing room temperature 
constituent (fiber, matrix, and interphase) properties of various materials. 

The primary purpose of the databank is to relieve the user of the repetitive 
burden of including this information as part of the primary input file. By 
assigning a unique code name for each constituent in the databank, the user 
merely has to specify the code name in the primary input file to acquire the 
numerous properties for a specific constituent. The assignment of unique 
code names also makes it convenient to update the databank as material 
properties become available. The databank is organized into three sections; 
the first for fiber properties, the second for matrix properties, and the 
third for interphase properties. Figure 1. 4.2-1 shows the databank organization 
and the properties contained for each constituent. Details on preparing the 
constituent databank are found in Section 3.0. 



Figure 1. 4.2-1 — Organization of the constituent databank. 
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1.4.3 


Primary Output File 


Results of the METCAN simulation are generated to a primary output file. 

The user has the flexibility to tailor the primary output file by specifying 
various combinations of the available output requests as part of the primary 
input file (see the PRINT input data record in Section 2.0). The simulation 
results comprise information at three levels, the laminate, ply, and 
constituent levels. The three levels of the analysis and the various 
information associated with each level are defined in Figure 1.4. 3-1. 



Figure 1.4. 3-1 — Levels of analysis information available in the primary output file. 
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1.4.4 


Postprocessing Output Files 


The ten postprocessing output files are optionally generated at the user’s 
request (see the POST input data record in Section 2.0). These ten files 
essentially contain the same information generated in the primary output 
file, but in a format from which selected output quantities can be extracted by 
the user and imported into graphics packages for convenient representation. The 
ten postprocessing output files and their contents are described in Figure 
1.4. 4-1. A more detailed description regarding the ten postprocessing files can 
be found in Section 4.0. 



Figure 1 .4,4-1 — Contents of the ten postprocessing files. 
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1 . 4.5 


Input/Output User’s Information 


The METCAN source code consists of approximately 7400 lines of FORTRAN 
statements. The program has been compiled using a variety of FORTRAN 77 
compilers on hardware ranging from microcomputers to supercomputers with 
relative ease. It is not anticipated that the user will have any difficulty 
with program installation. 

The file system for METCAN consists of: 


I) 

The primary input file 



FORTRAN 

I/O 

Unit 

5 

II) 

The constituent databank 

input file 

FORTRAN 

I/O 

Unit 

8 

III) 

The primary output file 



FORTRAN 

I/O 

Unit 

6 

IV) 

The ten postprocessing output 

files: 






i) 

Postprocessing 

File 

1 

FORTRAN 

I/O 

Unit 

31 


2) 

Postprocessing 

File 

2 

FORTRAN 

I/O 

Unit 

32 


3) 

Postprocessing 

File 

3 

FORTRAN 

I/O 

Unit 

33 


4) 

Postprocessing 

File 

4 

FORTRAN 

I/O 

Unit 

34 


5) 

Postprocessing 

File 

5 

FORTRAN 

I/O 

Unit 

35 


6) 

Postprocessing 

File 

6 

FORTRAN 

I/O 

Unit 

36 


7) 

Postprocessing 

File 

7 

FORTRAN 

I/O 

Unit 

37 


8) 

Postprocessing 

File 

8 

FORTRAN 

I/O 

Unit 

38 


9) 

Postprocessing 

File 

9 

FORTRAN 

I/O 

Unit 

39 


10) 

Postprocessing 

File 

10 

FORTRAN 

I/O 

Unit 

40 
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1.5 METCAN Flowchart 


The METCAN flowchart provides details of the execution sequence for a typical 
analysis and is shown in Figure 1.5-1. The analysis is initiated through the 
two input files. The primary input file contains the ply data, loading history, 
output requests, and specifies the constituent properties to be retrieved from 
the second input file, the constituent databank. Once the data from these two 
input files are read into METCAN, the code determines the ply thermal and 
mechanical properties followed by a laminate analysis to find the composite 
properties for the first increment. From this point, METCAN proceeds to the 
incremental iterative portion of the analysis. The microstresses in the 
various sub-regions of the unit cell are then calculated. These results are 
used in turn to update the constituent mechanical properties, the constituent 
thermal properties, and the constituent strength limits. The updated 
constituent properties are used in conjunction with micromechanics theory to 
determine the ply thermal and mechanical properties. Subsequently, a laminate 
analysis is performed to find the composite properties. At this point 
convergence is checked. If the convergence criteria is not met, the previously 
detailed steps are repeated until convergence is satisfied. 

Once convergence is achieved, results from the analysis are generated to 
the two output files according to the user specified requests. The first 
output file actually consists of a group of ten files referred to as the 
postprocessing files. These files are optionally generated by the user and 
basically contain the same information found in the other output file, the 
primary output file. The primary output file can be tailored by the user to 
contain various combinations of the results generated in the analysis. 

After the output requests have been generated, the incremental iterative 
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analysis is performed for the next load increment and these results 
generated to the output files. This cycle is continued until all the load 
increments have been analyzed. 



Figure 1.5-1— METCAN flowchart. 
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1.6 


METCAN Constitutive Relationship and Micromechanics Model 


Two features of the theoretical basis of METCAN are (1) the multifactor 
interaction relationship to functionally represent the nonlinear behavior of 
metal matrix composites at the constituent level and (2) the generic unit 
cell to derive the micromechanics equations. These two aspects are 
described in the following sections. 
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1.6.1 Nonlinear Constitutive Material Model 

METCAN treats material nonlinearity at the constituent level through the 
incorporation of a multifactor interaction relationship (MFIR). The MFIR 
models the material behavior in the form of a time-temperature-stress dependence 
of the constituent properties. The MFIR and the rationale for its utilization 
are expressed in Figure 1.6. 1-1. The form of the MFIR is a product function of 
terms raised to specified powers that define the values of the material property 
(P) of the fiber, matrix, and interphase. The constituent properties (P) which 
are determined through the MFIR are modulus, Poisson's ratio, strength, 
coefficient of thermal expansion, and heat capacity. For each fiber, matrix, 
and interphase property, there is generally a distinct set of values of the 
exponents necessary to express the functional definition of the property. The 
various exponents are determined from the appropriate experimental data. If 
there is a lack of experimental data, default values are used as an 
approximation. The values for the exponents are defined in the constituent 
databank for the various fibers, matrices, and interphases. The default values 
of the exponents are listed under the EXPONENTS record of each constituent in 
Section 3.0. 
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Rationale: 



• Gradual effects during most range, rapidly degrading near final stages 

• Representative of the in situ behavior for fiber, matrix, interphase, coating 

• Introduction of primitive variables (PV) 

• Consistent in situ representation of alt constituent properties In terms of PV 

• Room-temperature values for reference properties 

• Continuous interphase growth 

• Simultaneous interaction of ail primitive variables 

• Adaptability to new materials 

• Amenable to verification inclusive of alt properties 

• Readily adaptable to incremental computational simulation 



Notations: 

P - property; T - temperature; S - strength; R - metallurgical reaction; N - number of cycles; D - 
diameter; <r - stress; t - time; over dot - rate; subscripts: 0 - reference; F - final; M - mechanical; T - 
thermal 


Figure 1 .6.1-1 —Assumed multifactor Interaction relationship to represent the various factors which 
influence in Situ constituent materials behavior. 
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1 . 6.2 


Micromechanical Unit Cell and Subregions 


There are two alternative versions of the generic unit cell model used in 
METCAN, as shown in Figure 1.6. 2-1. A typical unit cell consists of a fiber and 
a matrix with or without an interphase. The unit cell is further subdivided 
into two or three subregions depending on the presence of an interphase. If an 
interphase is present, subregion A consists entirely of matrix material, 
subregion B consists of matrix and interphase, and subregion C consists of 
fiber, matrix, and interphase. If there is no interphase, two subregions exist. 
Subregion A consists entirely of matrix material, while subregion C consists of 
fiber and matrix, and subregion B does not exist. 


A: Matrix 

B: Matrix and interphase 3 

C: Matrix, fiber, and interphase 


Regions of 
constituent 
material 

non uniform tty — 


1 


Matrix 

' Interphase 



(a) Unit ceil with an interphase. 



(b) Unit ceil without an Interphase. 
Figure 1 .6.2-1 — METCAN generic unit cells. 
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1.7 


Notation and Units 


The notation used in the subsequent sections along with the corresponding 
units of measure are explained below. In principle, any consistent system 
of units may be used. However, the user should note that the sample input 

files and the sample constituent databank provided with METCAN are based on 

the given units. 


Symbols 

k 

t 

fc i 

z 

C 

C 

P 

D f 

E 

G 

K 

M 

N 

N f 

N M 


U 


Q 

S 

T 


Volume fraction 
Time 

Thickness 

Distance from reference plane 

Stress-strain relationship tensor 

Heat capacity 

Fiber diameter 

Normal modulus 

Shear modulus 

Ther ma 1 conduc t i v i ty 

In-plane moment resultant 

In-plane stress resultant 

Number of fibers per bundle 

Number of mechanical load cycles 

Number of thermal load cycles 

Lower surface pressure 

Upper surface pressure 

Out-of-plane shear resultant 

Strength 

Temperature 


sec 

in 

in 

psi 

Btu/lb 

in 

psi 

psi 

Btu/hr/in/°F 

lb-in/in 

lb/in 


psi 

psi 

lb/in 

psi 

o 

F 
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T 

a 


M 


MDE 

e 


v 


P 

<F 


'H 


0 


Superscripts 


Subscripts 


c 

d 

f 

i 

1 

m 

v 


x,y,z 
A, B, C 
T,C,S,x 
1 , 2,3 


Melting temperature 
Coefficient of thermal expansion 
Modified distortion energy 
Strain 

Poisson’ s ratio 

Density 

Stress 

Limit-state value of stress rate 
Orientation angle 


o 


F 


in/in/°F 


% 


lb/in 3 

psi 

psi/sec 

degrees 


Derivative with respect to time 


Composite (laminate) related quantity 

Interphase related quantity 

Fiber related quantity 

Index for i-th load/time step 

Lamina (ply) related quantity 

Matrix related quantity 

Void related quantity 

Global coordinate system axes 

Unit cell sub-regions 

Tension, compression, shear, and torsion, respectively 
Local coordinate system axes 
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2.0 Input Data Records 


The METCAN input file structure provides a straightforward user interface. The 
input file is organized into different records in a specific order. Thus each 
record must be included in the input file in the order they are defined in the 
following sections. Each record can be composed of one or several physical 
lines of data. Each line of data has a fixed format of ten eight-column fields 
(except for the title and comment records). The usual convention is that each 
record is identified by a character mnemonic in the first field of the data. 

The character mnemonics and other alphanumeric data are entered with character 
format (A8), while integer data are entered in integer format (18). Real data 
can be entered in either floating point (F8) or exponential (E8) formats. 
Character, integer, and exponential formats must be right justified, while 
the floating point format can be entered anywhere in the appropriate field. 
Figure 2.0-1 describes the individual records required in the input file, the 
mnemonic which identifies the record, the number of lines of data which comprise 
each record, and the order in which each record is read by METCAN. A 
description and example of each record is given in the following sections. 
Samples of complete input files are provided at the end of this section for 
illustration. 
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No mnemonic - 1 line 

Execution options 

1 . POST - 1 line 

2. LDDIST - 1 line 

Pty -description 

PLY - 1 line for each ply 

Composite material system 

MATCRD - 1 line for each material 

Cumulative load cycles 

CYCLES - 1 line 

Output report options 

1 . PRINT or OPTION - 1 or more lines 

2. PRINTOPT - 1 line 

Fiber-matrix Interphase 

INTRFACE - 1 line 

Analysis load history 

TMLOAD 

1 . Generalized table form - 3 or more lines 
or 

2. Simplified table form - 6 or more lines 


Figure 2.0-1 — Composition of the primary Input data file. 
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2.1 Input Data Record: $ 


Description : Comments 


Format and Example : 


$ COMMENT 

$ Any comment, up to 79 

characters 

long. 




Field No. 

Field Name 

Format 

Columns 

Description 



1 

$ 

A1 

i 

Mnemonic name 



2 

COMMENT 

A79 

2-80 

Comment 



Remarks : 








1. Comments can be inserted anywhere in the input data file. 
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2.2 Input Data Record: (No Mnemonic) 


Description : Problem title 

Format and Example : 

+ + + + + + h + h + 

TITLE 

+ + 1 + + + + + + + 

Any problem title, up to 80 characters long. 

+ + + + + + H v + + + 

Field No . Field Name Format Columns Description 
1 TITLE A80 1-80 Problem title 


Remarks : 

1. This record has no mnemonic. 

2. The first non-comment line of the input file is read as the title record. 
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2.3 Input Data Record: POST 


Description: Postprocessing files generation option 


Format and 

Example: 

-4- 






POST 

ID 

¥ 

4- 

T 







T 

T 



-i_ _ T - 

+ 

POST 

T 












-4-— » 

+ _ 




Field No. 

Field Name 

Format 

t 

Columns 

Description 



1 

POST 

A8 


1- 8 

Mnemonic name 



2 

ID 

A8 


9-16 

File generation 

option 



Remarks : 

1. There are two choices for the ID: 

T will result in the generation of ten postprocessing files 
F will not generate the postprocessing files 
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2.4 Input Data Record: LDDIST 


Description: Load redistribution option 

Format and Example : 


LDDIST 

ID 

- + + “ 



LDDIST 

F 




Field No. 

Field Name 

Format 

Columns 

Description 

1 

LDDIST 

A8 

1- 8 

Mnemonic name 

2 

ID 

A8 

9-16 

Load redistribution option 


Remarks : 

1. The load redistribution algorithm employs post failure models to simulate 
the progression of damage in the composite. Accumulated stresses are 

redistributed among the constituents at each load step. 

2. This option is useful when simulating composites containing highly 
nonlinear matrix materials. 

3. There are two choices for the ID: 

T will invoke the redistribution option 
F will not invoke the redistribution option. 
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2.5 Input Data Record: PLY 


Description: Defines each ply of the composite 


Format and Example : 


PLY 

PID 

— +. 
MID 

e i 


— + + + 

-+ 

— f — 

PLY 

1 

1 

45. 

0.005 




Field No. 

Field Name Format 

Columns 

Description 


1 

PLY 


A8 

1- 8 

Mnemonic name 


2 

PID 


18 

9-16 

Ply identification number 


3 

MID 


18 

17-24 

Material system identification 
number 

4 

e i 


F8 

25-32 

Orientation angle of the 

ply 

5 



F8 

33-40 

Ply thickness 


Remarks: 







1 . Each ply of the 

composite 

is defined by 

a separate PLY record. 



2. The ply ordering sequence is shown in Section 1.3.3. 

3. Each different MID must have a corresponding MID on a MATCRD record. 

4. The orientation angle is taken as the angle of the ply material 
coordinate system with respect to the structural coordinate system, 
shown in Section 1.3.4. 

5. If the tj field is left blank, the thickness is defaulted to 0.005 inches. 
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2.6 Input Data Record: MATCRD 


Description: Defines a composite material system (i.e. fiber/matrix 

composition) 

Format and Example : 


+ + + + + + + + + + 

MATCRD MID k_ k CODE 

f v 

+ + + + + + + + + + 


MATCRD 

1 0.35 

+ 

0.0 

+ 

SICATI 15 

Field No. 

Field Name 

Format 

Columns 

1 

MATCRD 

A8 

1- 8 

2 

MID 

18 

9-16 

3 

k f 

F8 

17-24 

4 

k 

V 

F8 

25-32 

5 

CODE 

A8 

33-40 


Description 
Mnemonic name 

Material system identification 
number 

Fiber volume fraction 
Void volume fraction 
Fiber and matrix code names 


Remarks : 

1. A separate MATCRD record must be defined for each different material system 
in the composite. 

2. Each different MID field must have a corresponding MID on a PLY record. 

3. The fiber volume fraction is restricted to the range: 

0.05 < k f < 0.78 

4. The CODE field specifies the four letter code names of the fiber and 
matrix and identifies the corresponding entries in the constituent 
databank. 
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2.7 Input Data Record: CYCLES 



CYCLES 

n m 

N I 

N N N N N N 

MFf TFf MFm TFm MFi TFi 

CYCLES 

1000. 

10. 

1.0E+06 4.0E+02 1.0E+06 4.0E+02 1 . 0E+06 4.0E+02 


Field No. Field Name' Format Columns Description 


1 

CYCLES 

A8 

1- 8 

Mnemonic name 

2 

n m 

F8 

9-16 

Number of mechanical cycles 

3 

n t 

F8 

17-24 

Number of thermal cycles 

4 

N MFf 

F8 

25-32 

Fiber mechanical cycles limit 

5 

N TFf 

F8 

33-40 

Fiber thermal cycles limit 

6 

N MFm 

F8 

41-48 

Matrix mechanical cycles limit 

7 

N 

TFm 

F8 

49-56 

Matrix thermal cycles limit 

8 

N 

MFi 

F8 

57-64 

Interface mechanical cycles limit 

9 

N Tr . 

TFi 

F8 

65-72 

Interface thermal cycles limit 

Remarks : 

1. The 

CYCLES record 

specifies 

an option 

to degrade the material properties 


due to cumulative mechanical and thermal load cycles. 

2. Specifying a value of 1.0 for N^ nullifies mechanical cyclic effects. 

Similarly, a value of 1.0 for N^ nullifies thermal cyclic effects. A value 
of 1.0 only nullifies the specified cycling. 
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3. 


The limit state values of the mechanical and thermal cycles (T^p and 

T ) can be specified for each constituent on this record. These fields 
TF 

are optional. 

4. The default limit state value of the mechanical cycle (N^-J is set 
at 10 6 cycles. 

5. The default limit state value of the thermal cycle (N^.p) is set internally 
at 400 cycles. 

6. In the case where both cumulative mechanical and thermal cycles and a 
linear load history are desired additional fields are required on the 
TMLOAD (Simplified Table Form) record to differentiate between cyclic and 
static loads. See the TMLOAD (Simplified Table Form) record. 

7. When cycling is specified, output information for the cycling portion of 
the load history is generated only for the last cycle. 
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2.8 Input Data Record: PRINT 


Description: Defines output request options 

Format and Example : 

+ + + + + + + + + + 

PRINT OUTCODE 

+ + + + + +- + + + + 

PRINT MICRO 

+ + + h + + + + + + 

Alternate Format and Example : 


+_ 

OPTION 

OUTNUM 








OPTION 
+ - 

10 









Field No. 

Field Name 

Format 

Columns 

Description 

1 

PRINT 

A8 

1- 8 

Mnemonic name 

2 

OCODE 

A8 

9-16 

Output option request code 


Alternate Form: 


Field No. 

Field Name 

Format 

Columns 

Description 

i 

OPTION 

A8 

1- 8 

Mnemonic name 

2 

OUTNUM 

18 

9-16 

Output option request number 
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Remarks: 


1. This record allows the user to tailor the output generation. 

2. There are two alternative forms to specify the output requests: 

(i) An output code form (OUTCODE) 

(ii) An output number form (OUTNUM) 

The output code form is recommended, since the output number form will be 
phased out in subsequent versions. 

3. The available output request options are: 


Output Code 
OUTCODE 

Output No. 
OUTNUM 

Description 

ALL 

0 

Complete report of all available information 

CONST I 

8 

Laminate constitutive relationships 

DISPFOR 

11 

Displacement-force relationships 

FEMDATA 

10 

Information for finite element analysis 

FLINDEX 

5 

Report of failure index 

LDSTEP 

29 

Information for current load/time step 

MICRO 

28 

Constituent total stresses and strains for each 
ply 

NONE 

99 

Output reporting suppressed 

PLYRESP 

12 

Ply properties and response variables 

PLYSTRS 

15 

Ply stresses and strains 

PROPCOM 

7 

Laminate 2-D and 3-D properties 

PROPCUR 

27 

Constituent instantaneous properties 

PROPREF 

26 

Constituent reference properties 

REDSTIF 

9 

Laminate reduced membrane and bending stiffnesses 
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STRCON 16 

STRSINF 18 

STRSTRN 6 


Stress concentration factors around a circular 
hole in an infinite plate 
Stress-strain influence coefficients 
Laminate stress-strain relationships and 
MSC/NASTRAN MAT9 record 


4. There is no limit to the number of output options which can be requested, 
however each output option requires a separate record. 

5. The output reporting frequency is controlled through the PRINTOPT record. 
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2.9 Input Data Record: PRINTOPT 


Description: Defines the output reporting frequency 

Format and Example : 

+ + + + + + + + + + 

PRINTOPT PR ID 

PRINTOPT 10 

Field No . Field Name Format Columns Description 

1 PRINTOPT A8 1- 8 Mnemonic name 

2 PRID A8 or 18 9-16 Output reporting frequency 

Remarks : 

1. This record is optional. The default is to print the output for each load 
step of the analysis. 

2. There are three choices for the PRID: 

(i) ALL Output is produced for each load step of the analysis 

(ii) LAST Output is produced for only the last load step of the 

analysis 

(iii) i Output is produced after every i-th step of the analysis 
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2.10 Input Data Record: INTRFACE 


Description: Defines an interphase between the fiber and matrix 

Format and Example : 


+— 

+-- 


,+ 

-+- 

+ 

+ + + 

+ 

INTRFACE 

INID 

4 d 






+ 








INTRFACE 

1 

0.10 














Field No. 

Field 

Name 

Format 


Columns 

Description 


1 

INTRFACE 

A8 


1- 8 

Mnemonic name 


2 

INID 


18 


9-16 

Interphase indicator 


3 



F8 


17-24 

Interphase thickness 


Remarks: 









1. There are two choices for the INID field: 

(i) 0 An interphase is not considered in the analysis 

(ii) 1 An interphase is considered in the analysis 

2. If an interphase is not considered (INID=0), then perfect bonding between 
the fiber and matrix is assumed. 

3. The field t^ is only valid for an INID choice of 1 and specifies the 

thickness of the interphase as a fraction of the fiber diameter. The 

thickness performs the same function as the fiber and matrix code 
names in the MATCRD record. It specifies the interphase entry in the 
constituent databank containing the room temperature material properties 
of the interphase. 

4. Only one interphase can be specified in the input data file. 
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2.11 Input Data Record: TMLOAD (Generalized Table Form) 


Description: Defines analysis load history with discrete points 

Format and Example : 


+ — 

TMLOAD 

NLS 








TMLOAD 

4 









*1 T 1 

1,1 

T 

1 2, i 


.... 














0.0 2000. 

1800. 


.... 





. 

+ . 

+ 

+ 

+ 

+ 

+ 

+ 

+ 



T 

1 k,i 


.... 



T 

1 n, i 

+ 

+ +_ 

+ 

_ + 





1200. 


.... 



1000. 


+ + + + + + + + + + 


N , 
x, 1 

N , 
y. i 

N , 
xy, i 

M , 
x, i 

M . M , Q . Q , P. . P. 

y,l xy.l xz,i yz,i 1,1 u,i 

100. 

100. 

0 . 

150. 

150. 0. 0. 0. 2000. 2000. 


+ + + + + + + + + + 
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Field No. 

Field Name 

Format 

Columns 

Description 

Line 1: 

1 

TMLOAD 

A8 

1- 8 

Mnemonic name 


2 

NLS 

18 

9-16 

Number of load history points 

Line 2: 

1 

t. 

X 

F8 

1- 8 

Value of time 

2 

T 

1 1 * i 

F8 

9-16 

Temperature of ply 1 at time t^ 

3 

T 

1 2 , i 

F8 

17-24 

Temperature of ply 2 at time t^ 

4-10 

T -T 

1 3, i 1 9 , i 

6(F8) 

25-80 

Temperature of plies 3 through 9 at 
time t. 

l 

Line 3: 

1-10 

T i 

1 n, i 

10 (F8 ) 

1-80 

Temperature of plies 10 through n 
at time t^ (see Remark 4) 

Line 4: 

1-3 

N (x,y,xy), i 

3(F8) 

1-24 

Global in-plane force resultant 
at time t. 

l 

4-6 

M (x,y,xy) , i 

3(F8 ) 

25-48 

Global in-plane moment resultant at 
time t^ 

7-8 

^(xz.yz) , i 

2(F8) 

49-64 

Global transverse shear resultant 
at time t. 

l 

9-10 

P ( 1 , u) , i 

2 (F8 ) 

65-80 

Global lower and upper surface 
pressure at time t^ 
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Remarks: 


1. There are two alternative forms of the TMLOAD record, which are specified 
by the sign of NLS. 

(i) A positive NLS defines the loading history as a number of 
discrete points, record TMLOAD (Generalized Table Form). 

(ii) A negative NLS defines a number of segments to define a piecewise 
linear loading history, record TMLOAD (Simplified Table Form). 

2. Line one, containing the mnemonic TMLOAD, is specified only once in the 
input data file. 

3. Lines two onward are repeated for every discrete point specified (NLS 
times) . 

4. The second line of data specifies the time t^ and the temperature of each 
ply at this time. If there are more than nine plies, the temperatures 
for the tenth ply onwards can be continued on subsequent lines until all 
plies are represented. 

5. Ply temperatures are restricted to values below the lowest melting 
temperature of the constituents specified on the MATCRD and INTRFACE 
records. 

6. The last line of data specifies the various mechanical loads applied on 
the composite at the given instant. The positive sign conventions for 
the mechanical loads are shown in Section 1.2.5. 

7. Default values for the temperatures and mechanical loads are zero. 

8. An example of this record, along with a figure, can be found in Case 2 of 
the sample input data files (Sec. 2.13). 
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2.12 Input Data Record: TMLOAD (Simplified Table Form) 

Description: Defines analysis load history with linear segments 


Format and Example: 


TMLOAD NLS CYCPRO 



t t i+1 INCH CYCOPT 

+ + + + 


0.0 300. 35 BEGCYCL UPPER 



T 

1 1,1+1 

T 1 k, i+1 

T 1 n, i+1 

70. 

70. 

70. 
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.+ + + + + + + + + + 


N « 
x,i 

N < 
y. i 

N < 
xy, i 

M i 

x, i 

V Vi Q xz.i Vi 

p i.i 

P u, i 

100. 

100. 

0 . 

150. 

150. 0. 0. 0. 

2000. 

2000. 


.+ + + + + + + + + + 


+ + + + + + + + + + 

^x,i+l ^y, i+i ^xy,i+l ^x,i+l ^y, i+1 ^xy,i+l ^xz,i+l ^yz,i+l P l»i+1 P u,i+1 
+ + + + + + + + + + 

0 . 0 . 0 . 0 . 0 . 0 . 0 . 0 . 0 . 0 . 

+ + + + + + + + + + 


Field No. 

Field Name 

Format 

Columns 

Line 1: 




1 

TMLOAD 

A8 

1- 8 

2 

NLS 

18 

9-16 

3 



(skip field 

4 

CYCPR0 

A8 

25-32 


Description 
Mnemonic name 

Number of load history segments 

Alphanumeric to indicate processing 
when cyclic loads are requested 
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Line 2: 


1 

h 

F8 

00 

1 

H 

Value of time at t^ 

2 

h+i 

F8 

9-16 

Value of time at t^ + ^ 

3 

I NCR 

18 

17-24 

Number of increments between time 


CYCOPT 


h and t i + l 


(skip field 4) 


33-48 Alphanumeric to differentiate 

between cyclic and static loadings 


1 T M * 

1 1,1 

T F8 

1 2, i 

T, 0 . -T. . 7 (F8) 

13,1 1 n, l 


1 - 8 
9-16 
17-80 


Temperature of ply 1 at time t^ 
Temperature of ply 2 at time t^ 
Temperature of plies 3 through n at 


time t. 


1 

T i 

l.i+i 

F8 

1- 8 

Temperature 

of 

ply 1 

at time t. , 
l + l 

2 

T i 

2, i+1 

F8 

9-16 

Temperature 

of 

ply 2 

at time t. _ 
l+l 

3-10 

T i 

3, i+1 

7 (F8) 

17-80 

Temperature 

of 

plies 

3 through n at 


1 n, i + 1 


time t. 


N ( . . 3(F8) 

(x,y,xy) , l 


M, , . 3(F8) 

(x,y,xy) , i 


1-24 Global in-plane force resultant at 
time t^ 

25-48 Global in-plane moment resultant at 
time t. 
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7-8 


2 (F8 ) 


(xz,yz) , i 


9-10 P ri , . 2 (F8 ) 

(1 ,u ) , 1 


49-64 Global transverse shear resultant 

at time t^ 

65-80 Global lower and upper surface 

pressure at time t^ 


Line 6: 





1-3 

N ( x , y , xy ) , i+1 

3(F8) 

1-24 

Global in-plane force resultant 
at time t^ + ^ 

4-6 

M (x,y,xy) , i+1 

3(F8) 

25-48 

Global in-plane moment resultant 

time t. . 
l+l 

7-8 

^(xz.yz) , i+1 

2(F8) 

49-64 

Global transverse shear resultant 
at time t^ + ^ 

9-10 

P ( 1 , u) , i+1 

2 (F8 ) 

65-80 

Global lower and upper surface 
pressure at time t^ + ^ 


Remarks: 


1. There are two alternative forms of the TMLOAD record, which are specified 
by the sign of NLS. 

(i) A positive NLS defines the loading history as a number of 
discrete points, record TMLOAD (Generalized Table Form). 

(ii) A negative NLS defines a number of segments to define a piecewise 
linear loading history, record TMLOAD (Simplified Table Form). 

2. The maximum value of NLS is -10 for the Simplified Table Form. 
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3. The field CYCPRO on the first line of data is only required if cumulative 
mechanical and thermal cycles are specified on the CYCLES record. There 
are two choices for this field: 

(i) PROCESS Indicates that processing is simulated in the 

loading history 

(ii) NOPROCES Indicates that processing is not considered. 

This is the default choice. 

4. Lines two onward are repeated NLS times. 

5. The second line of data specifies the times and t^ + ^ f the number 
of linear time segments into which the two times are divided, and 
differentiates between cyclic and static loadings. 

6. The field CYCOPT in the second line of data is only required if cumulative 
mechanical and thermal cycles are specified on the CYCLES record. There 
are five choices for this field: 


(i) 

BEGPROCE 

Denotes 

the 

beginning of processing 

(ii) 

ENDPROCE 

Denotes 

the 

end of processing 

(iii) 

BEGCYCL UPPER 

Denotes 

the 

initial portion of cycling 

(iv) 

ENDCYCL LOWER 

Denotes 

the 

final portion of cycling 

(v) 

ENDCYCL 

Denotes 

the 

end of cycling 


These alphanumerics are required to differentiate between the cyclic and 
static loads. 

7. Ply temperatures are restricted to values below the lowest melting 

temperature of the constituents specified on the MATCRD and INTRFACE 
records . 
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8. As a general rule in choosing the number of increments (INCR), select INCH 
such that each subincrement represents a change in temperature of about 
35°F. For example, if T. . - 1470 F and T = 70 F, then INCR = 

i + l i 

( 1470-70 )/35 = 40 increments. If this INCR turns out to be less than 20, 
set the INCR to 20. 

9. The third line of data specifies the ply temperatures at time t^. If there 
are more than ten plies, the temperatures can be continued on as many 
subsequent lines as necessary. Once all the ply temperatures for time 

t^ are given, the ply temperatures at time t^ + ^ are specified on the next 
line of input, in a similar fashion. 

10. The last two lines of data specifies the various mechanical loads applied 
on the composite at times t^ and t. + ^, respectively. The positive sign 
conventions for the mechanical loads are shown in Section 1.1.5. 

11. Default values for the temperatures and mechanical loads are zero. 

12. A static example of this record can be found in Case 1 of Section 2.13, 
while a cyclic example can be found in Case 2 of Section 2.13. 
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2.13 Sample Input Data Files 

Three complete input data files, highlighting different features of METCAN are 
provided below for illustration. The first sample input data file is a more 
detailed version of the sample problem described in Section 1.2, with each data 
record discussed in depth. Comment cards are also included throughout the first 
sample input data file to explain the function of each data record. The second 
and third sample input data files point out different features of METCAN, but 
are not presented in the same detail as the first sample input file. 
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The first sample input data file (Figure 2.13-1) simulates a two ply 
unidirectional SiC/Ti-15-3 composite under a longitudinal tensile load at room 
temperature. The first line of the sample input file is the title card. This 
record contains no mnemonic and can contain any title up to 80 characters long. 
The second line of the sample input data file is a comment card, denoted by ’ $* 
in the first column. Comment cards are inserted before each different record of 
the sample input data file to briefly describe the function of each record. The 
next record, containing the mnemonic POST, is used to control the generation of 
the ten postprocessing files. In the sample input, the 'F* in the second field 
denotes that these files are not to be generated. The next record (containing 
the LDDIST mnemonic) controls the load redistribution option in METCAN. Similar 
to the POST record, the * F f in the second field specifies that this option is 
not to be employed. A * T* in the second field of either the POST or LDDIST 
records would utilize these options. 

The next record begins with the mnemonic PLY. This record consists of five 
fields and contains information regarding the ply configuration. Fields two 
through five contain the ply number, the corresponding material system number 
used in the MATCRD record, the orientation angle of the ply, and the ply 
thickness, respectively. A separate PLY record is required for each ply in the 
laminate. For the case depicted in the sample input data file, two PLY records 
are required since the laminate is composed of two plies. Each ply has an 

o 

orientation angle of 0 , a thickness of 0.005 inches, and contains the same 
materials (represented by the same material system number). 

The MATCRD record contains four fields and includes the material details. The 
second field contains the material system number corresponding to the third 
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SAMPLE INPUT DATA FOR METCAN, CASE 1 
S No postprocessing files requested. 

POST F 

$ Load redistribution option turned off. 

LDDIST F 

. S Ply configuration details. 

PLY 1 1 0.0 0.005 

PLY 2 1 o.O 0.005 

S Material details. 

HA" 1 ™ 0 1 .30 0. S I CAT 115 

S Mechanical and thermal cycling effects. 

CYCLES 1.0 l.o 

$ Output requests. 

PRINT ALL 

$ Print output for every load step. 

PRINTOPT ALL 

$ Interface details. 

INTRFACE 0 . io 

$ Simplified table input. 

TMLOAD -2 

$ Start time, end time, and number of increments. 

0.0 1000. 35 

S Temperature in each ply at the beginning and end of the time step. 
1200. 1200. 

70. 70. 

S Mechanical loads at the beginning and end of the time step. 

°’ °* °- 0 . 0 . ' 0 . 0 . 0 . 0 . 
0> °* °* 0* 0 . 0 . 0 . o. 

t Start time, end time, and number of increments. 

1000. 2000. 25 

? Tem P erature in each ply at the beginning and end of the time step. 
70. 70. 

70. 70. 

Mechanical loads at the beginning and end of the time step. 

°' °- °’ 0 . o. 0 . 0 . o. 

0- 0. 0. 


500. 0. 0. 0. 

+ ♦ + +_ 


0 . 0 . 

+ ♦ 


Figure 2.1 3-1 —Sample case 1 Input data file. 
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field of the PLY record. Each existing material system number in the PLY 
record must have a matching MATCRD record with the same material system 
number. Fields three through five respectively contain the fiber volume ratio, 
the void volume ratio, and the fiber and matrix codes used to identify the 
corresponding entries in the constituent databank. For the sample input data 
file, this record shows that each ply contains a fiber volume ratio of 0.30, a 
void volume ratio of 0.0, and is composed of a silicon carbide fiber and a 
titanium (Ti-15-3) matrix. 

The next record begins with the mnemonic CYCLES and is comprised of data 
regarding thermal and mechanical cycling. Fields two and three specify the 
number of mechanical and thermal cycles that the laminate is subjected to. A 
value of 1.0 in either one of these fields nullifies the specified cycling. 

The sample input data file specifies no mechanical or thermal cycling. 

The PRINT and PRINTOPT records are used to tailor the output generation. The 
'ALL* in the second field of the PRINT record results in the generation of all 
the available output. The ’ALL’ in the PRINTOPT record defines that the 
output chosen in the PRINT record is printed at every load/time step of the 
analysis. 

The INTRFACE record is used to define an interface. The value in the second 
field determines whether an interface is present between the fiber and matrix or 
a perfect bond exists between the fiber and matrix. A '0* indicates that there 
is no interface present, while a * 1* denotes the presence of an interface. The 
third field is only used when an interface is present and specifies the 
thickness of the interface as a fraction of the fiber diameter. This value is 
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used to identify the interface entry in the constituent databank. For the 
sample input data file, there is no interface present. 

The next record, containing the mnemonic TMLOAD, defines the loading history. 

The second field, which contains ’-2’, specifies a linear loading history 
divided into two segments. The loading history is shown in Figure 2.13-2. Each 
segment of the loading history is represented by five lines of data which contain 
time, temperature, and loading conditions. For the first segment, the first 
line of data consists of three fields which contain the initial and final times 
for the segment and the number of increments into which the segment is divided. 
The second line consists of two fields and contains the temperature of each ply 
at the initial time of the segment, while the third line contains the 
temperature of each ply at the final segement time. The mechanical loads are 
specified in lines four and five for the initial and final segment time in a 
fashion similar to the temperatures. The second segment of the loading history 
is defined in a similar manner with the same group of five lines. In the sample 
input data file, the first segment extends from time 0 to 1000 seconds and is 
divided into 35 increments. Both plies are cooled from 1200°F to 70°F in the 
absence of any mechanical loads. The second segment continues from 1000 to 2000 
seconds and is divided into 25 increments. Each ply is held at a constant 
temperature of 70 F and a 500 lb/in force resultant is applied in the global 
x-direction. 
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Mechanical load, Ib/ln Temperature, 




Figure 2.13-2 — Sample case 1 loading history: simplified table 
toon for static loading. 
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The second sample input data file (Figure 2.13-3) models a two ply [90/0] W/Cu 
composite with a total thickness of 0.20 inches and a fiber volume ratio of 
0.35. Postprocessing files are requested and the load redistribution option is 
not in effect. Mechanical and thermal cycling is turned off. Three output 
options are requested, 1) information regarding the current load/time step, 2) 
laminate stress-strain relationships, and 3) constituent microstresses and 
microstrains. The output is reported after every fifth time step. An interface 
with a thickness of 2% of the fiber diameter is specified. The generalized 
table form is used to define the loading history in two discrete points as shown 
in Figure 2.13-4. The first point is defined at time 0 seconds, with the 
temperature of both plies set to 1500 F and a stress resultant of 1000 lb/in 
applied in the global x-direction. The second point occurs at time 500 seconds, 
the temperature of the two plies is reduced to 700 F and the stress resultant in 
the x-direction is decreased to 800 lb/in. For simplicity, only two points are 
defined in this example. Generally for the loading history shown, many more 
points would be required to accurately model the loading history. 
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SAMPLE INPUT DATA FOR METCAN, CASE 2 


POST 

T 



LDDIST 

F 



PLY 

1 

1 

90. 0.010 

PLY 

2 

1 

0. 0.010 

MATCRD 

1 

.35 

0. TUNGCOPR 

CYCLES 

1.0 

1.0 1. 

. 0E+O6 4.0E+02 1.0E+06 4.0E+02 1 . 0E+06 4.0E+02 

PRINT 

LDSTEP 



PRINT 

STRSTRN 



PRINT 

MICRO 



PRINTOPT 

5 



INTRFACE 

1 

.02 


TMLOAD 

2 



0.0 

1500. 

1500. 


1000. 

0. 

0. 

o 

o 

o 

o 

o 

o 

500. 

700. 

700. 


500. 

+- 

0. 

0. 

0. 0. 0. 0. 0. 0. 


Figure 2.1 3-3 — Sample case 2 input data file. 
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Figure 2.1 3-4— Sample case 2 la 
form for static loading. 
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The third sample input data file (Figure 2.13-5) corresponds to the first sample 
case with the addition of mechanical and thermal cycling as shown in Figure 
2.13-6. The addition of cycling is reflected in the CYCLES and TMLOAD records. 
The CYCLES record specifies 750,000 mechanical load cycles and 300 thermal load 
cycles. The TMLOAD record defines a piecewise linear loading history divided 
into four segments and indicates that processing will be simulated in the 
loading history. The first segment ranges from time 0 to 1000 seconds and is 
divided into 35 time steps. The alphanumeric BEGPR0CE denotes that the first 
segment simulates the processing of the composite by cooling down from 1200°F to 

o 

70 F in the absence of mechanical loads. The second segment ranges from 1000 to 
2000 seconds and is composed of 25 time steps. The BEGCYCL UPPER alphanumeric 
denotes that this segment is the initial portion of the cyclic loading, which is 
composed of a heat up from 70°F to 1000°F with the addition of a 250 lb/in force 
resultant in the global x-direction. The third segment spans time 2000 to 3000 
seconds and is split into 25 time steps. The alphanumeric ENDCYCL LOWER 
indicates that this segment is the final portion of the cyclic loading, which 
consists of a cool down from 1000°F to 70°F and an unloading of the force 
resultant in the global x-direction from 250 lb/in to zero. The fourth segment 
begins at 3000 seconds and ends at 4000 seconds. It is divided into 25 time 
steps and the alphanumeric ENDCYCL indicates that this segment is not to be 
cycled. The loading for this segment consists the application of a 500 lb/in 
force resultant in the global x-direction at a constant temperature of 70°F. 
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SAMPLE INPUT DATA FOR METCAN, CASE 3 


POST F 

LDDIST F 


PLY 

1 

1 

0.0 0.005 




PLY 

2 

1 

0.0 0.005 




MAT CRD 

1 

.30 

0. SICATI15 




CYCLES 

7.SE+05 3. 

. 0E+02 1. 

. 0E+06 4.0E+02 1.0E+06 

4.0E+02 1. 

. 0E+06 

4.0E+02 

PRINT 

ALL 






PRINTOPT 

ALL 






INTRFACE 

0 

. 10 





TMLOAD 

-4 

PROCESS 




0.0 

1000. 

35 

BEGPROCE 




1200. 

1200. 






70. 

70. 






0 . 

0 . 

0 . 

o 

o 

o 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

o 

o 

o 

0 . 

0 . 

0 . 

1000. 

2000. 

25 

BEGCYCL UPPER 




70. 

70. 






1000. 

1000. 






0 . 

0 . 

0 . 

o 

o 

o 

0 . 

0 . 

0 . 

250. 

0 . 

0 . 

o 

o 

o 

0 . 

0 . 

0 . 

2000. 

3000. 

25 

ENDCYCL LOWER 




1000. 

1000. 






70. 

70. 






250. 

0 . 

0 . 

o 

o 

o 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

o 

o 

o 

0 . 

0 . 

0 . 

3000. 

4000. 

25 

ENDCYCL 




70. 

70. 






70. 

70. 






0 . 

0 . 

0 . 

o’ 

o 

o' 

0 . 

0 . 

0 . 

500. 

0 . 

0 . 

o 

o 

o 

0 . 

0 . 

0 . 


+ + + ♦ + + — 


Figure 2.1 3-5 — Sample case 3 input data file. 
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3 . 0 Const i tuen t Databank 


The databank is organized into three sections with each section containing a set 
number of records in a predefined order. The first section contains the fiber 
properties and contains nine records. The second section arranges the matrix 
properties into ten records, while the third section contains the interphase 
properties on nine records. Each record is composed of one or more lines of 
data with a varying format from record to record. Note that the format for the 
constituent databank records is different from the fixed format employed in the 
primary input file. Each line of data begins with a mnemonic. Character, 
integer, and exponential formats must be right justified, while the floating 
point format can fall anywhere in the appropriate field. The various 
constituent properties input into the databank are room temperature values. A 
sample constituent databank is presented at the end of the section for 
illustration. 
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3.1.1 


Fiber Data Record: CODE 


Description : Defines fiber code name 


Format and Example : 


CODE 

COMMENT 









SICA 

Optional description of fiber 







Field No. 

Field Name 

Format 

Columns 

Description 

1 

CODE 

A4 

1-4 

Fiber code name 

2 

COMMENT 

A76 

5-80 

Optional description of fiber 

Remarks : 





1. The 

CODE for each 

fiber in 

the databank 

must be unique. 


3-4 


3.1.2 Fiber Data Record: FP 


Description : Defines fiber physical properties 


Format 

and 

Example: 





+-- 





FP 

N f 

D f 

p f 

T 

Mf 


_+ 

+— 





FP 

1 

. 560E-02 . 

110E+00 

. 487E+04 








Field 

No. 

Field Name 

Format 

Columns 

Description 

1 


FP 

A2 

1-2 

Mnemonic name 

2 


N f 

15 

4-9 

Number of fibers per bundle 

3 


D f 

E10.3 

10-19 

Fiber diameter 

4 


p f 

E10.3 

20-29 

Fiber weight density 

5 


T 

Mf 

E10.3 

30-39 

Fiber melting temperature 
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3.1-3 Fiber Data Record: FE 


Description : Defines fiber elastic properties 


Format and Example : 


+ — 

E fll 

+ 

F 

1 22 

+ — 

V fl2 

+ 

f23 

G f 12 

+ 

G f 23 

: . 620E+08 

. 620E+08 


■ V " 

. 300E+00 
+ 

. 300E+00 
+ 

. 238E+08 
+ — 

. 238E+08 


Field No. 

Field Name 

Format 

Columns 

Description 

1 

FE 

A2 

1-2 

Mnemonic name 

2 

E fll 

E10.3 

4-13 

Fiber 

longitudinal modulus 

3 

E f22 

E10.3 

14-23 

Fiber 

transverse modulus 

4 

V f 12 

E10.3 

24-33 

Fiber 

Poisson’ s ratio 

5 

V f23 

E10.3 

34-43 

Fiber 

Poisson’s ratio 

6 

G f 12 

E10.3 

44-53 

Fiber 

in-plane shear modulus 

7 

G f 23 

E10.3 

54-63 

Fiber 

out-of-plane shear modulus 
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3.1.4 Fiber Data Record: FT 


Descript ion : Defines fiber thermal properties 


Format and Example : 


FT 

-+- 


a fn 

+ . 

a f22 

K m 

+ _ 

K f22 

C pf 

FT 


, 180E-05 . 180E-05 

. 750E+00 

. 750E+00 

. 290E+00 




_____ + 





Field 

No. Field 

Name 

Format 

Columns 

Description 


1 

FT 


A2 

1-2 

Mnemonic name 


2 

a fn 


E10.3 

4-13 

Fiber longitudinal coefficient of 







thermal expansion 


3 

“f22 


E10.3 

14-23 

Fiber transverse coefficient of 







thermal expansion 


4 

K fll 


E10.3 

24-33 

Fiber longitudinal thermal 







conductivity 


5 

K f22 


E10.3 

34-43 

Fiber transverse thermal 







conductivity 


6 

C . 
Pt 


E10. 3 

44-53 

Fiber heat capacity 
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3.1.5 Fiber Data Record: FS 


Description: Defines fiber strength allowables 


Format and Example: 


FS S f 1 IT S fllC S f22T S f22C S fl2 S f23 

- + + + + + + + 

FS . 500E+06 . 650E+06 . 500E+06 .650E+06 . 300E+06 .300E+06 


Field No. 

Field Name 

Format 

Columns 

Description 

1 

FS 

A2 

1-2 

Mnemonic name 

2 

S f 1 IT 

E10.3 

4-13 

Fiber longitudinal tensile strength 

3 

S fllC 

E10.3 

14-23 

Fiber longitudinal compressive 





strength 

4 

S f22T 

E10.3 

24-33 

Fiber transverse tensile strength 

5 

S f22C 

E10.3 

34-43 

Fiber transverse compressive 





strength 

6 

s 

fl2 

E10.3 

44-53 

Fiber in-plane shear strength 

7 

s 

f 23 

E10.3 

54-63 

Fiber out-of-plane shear strength 


Remarks : 

1. In cases where only the experimental tensile strength is available, the 
following approximate values are used for the compressive and shear 
strengths: 

= S.j. and Sg = 0.6 S^, 
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3.1.6 Fiber Data Record: SIGFO 


Description : Defines reference values of fiber stress 


Format and Example : 


SIGFO 

°fll 

+ _ 

°f22 

-+ 

<r fl2 

0* 

f23 

+ - 

0 * 

f 13 

+ 

°f33 

SIGFO 

0 . 

+ - 

0 . 

0 . 

- + 

0 . 

+ - 

0 . 

+ — 

0 . 




Field No. 

Field Name 

Format 

Columns 

Description 

1 

SIGFO 

A5 

1-5 

Mnemonic name 

2 

11 

E10.3 

11-20 

Reference longitudinal stress value 

3 

°f22 

E10.3 

21-30 

Reference transverse stress value 

4 

cr 

f 12 

E10.3 

31-40 

Reference in-plane shear stress 





value 

5 

O' 

f23 

E10.3 

41-50 

Reference out-of-plane shear stress 





value 

6 

<r 

f 13 

E10.3 

51-60 

Reference out-of-plane shear stress 





value 

7 

cr 

f33 

E10.3 

61-70 

Reference through the thickness 


stress value 
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Remarks: 


1. The reference values of stress are used in the stress term of the 

multifactor interaction relationship and are typically set to zero for 
convenience. 
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3-11 


fiber stress rate 


O' & & 

f 23 f 13 f33 

+ + 

0 * 0 . 0 . 


Description 


Mnemonic name 

Reference longitudinal stress rate 
value 

Reference transverse stress rate 
value 

Reference in-plane shear stress 
rate value 

Reference out-of-plane shear stress 
rate value 

Reference out-of-plane shear stress 
rate value 

Reference through the thickness 


stress rate value 



Remarks: 


1. The reference values of stress rates are used in the stress rate term of 

the multifactor interaction relationship and are typically set to zero for 
convenience. 
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3 . 1.8 


Fiber Data Record: EXPONENTS 


Description : Defines values of exponents in the multifactor interaction 

relationship 

Format and Example : 


+. 

EXPONENTS 
+ - 

+. 

FEXl 

FEX2 

+ - 

FEX3 

FEX4 

FEX5 

FEX6 

+- 

FEX7 

FEX8 

FEX9 

FEX10 

EXPONENTS 

0.25 

0.25 

0.25 

0.25 

0.25 

- 

0.00 

0.25 

0.25 

0.25 

0.25 

. ' ■ ' » “ T — 

EXPONENTS 
+- 

FEX11 

FEXl 2 
+ - 

FEXl 3 
+ - 

FEXl 4 

FEXl 5 
+- 

FEXl 6 

T" 

FEXl 7 
+ - 

FEXl 8 

f - 

FEXl 9 
+ - 

FEX20 


EXPONENTS 0.25 0.25 0.25 0.25 0.00 0.25 0.25 0.00 0.25 0.25 

+ + + + + + + + + + +. 


+ . 

EXPONENTS 
+. 

FEX21 

+ - 

FEX22 
+_ 

FEX23 

+ - 

FEX24 

FEX25 

+ - 

FEX26 

FEX27 

FEX28 

+ - 

FEX29 

FEX30 

EXPONENTS 

0.25 
+ - 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 


EXPONENTS 

FEX31 














+ 


+ 

+ 

+ - 

" + 

+- 


EXPONENTS 

0.25 
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Field No. 

Field Name 

Format 

Columns 

Description 

Line 1: 

1 

EXPONENTS 

A9 

1-9 

Mnemonic name 

2 

FEX1 

F5.2 

11-15 

Reserved for future use 

3 

FEX2 

F5.2 

17-21 

Modulus-stress rate exponent 

4 

FEX3 

F5.2 

23-27 

Modulus-stress exponent 

5 

FEX4 

F5.2 

29-33 

Modulus-temperature exponent 

6 

FEX5 

F5.2 

35-39 

Strength-stress rate exponent 

7 

FEX6 

F5.2 

41-45 

Strength-stress exponent 

8 

FEX7 

F5.2 

47-51 

Strength-temperature exponent 

9 

FEX8 

F5.2 

53-57 

Poisson’s ratio-stress rate 
exponent 

10 

FEX9 

F5.2 

59-63 

Poisson’s ratio-stress exponent 

11 

FEX10 

F5.2 

65-69 

Poisson's ratio-temperature 
exponent 

Line 2: 

1 

EXPONENTS 

A9 

1-9 

Mnemonic name 

2 

FEX11 

F5.2 

11-15 

Reserved for future use 

3 

FEX12 

F5.2 

17-21 

Reserved for future use 

4 

FEX13 

F5.2 

23-27 

Reserved for future use 

5 

FEX14 

F5.2 

29-33 

Heat conductivity-stress rate 
exponent 

6 

FEX15 

F5.2 

35-39 

Heat conductivity-stress exponent 

7 

FEX16 

F5.2 

41-45 

Heat conductivity-temperature 
exponent 

8 

FEX17 

F5.2 

47-51 

CTE-stress rate exponent 
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9 

FEX18 

F5.2 

53-57 

CTE-stress rate exponent 

10 

FEX19 

F5.2 

59-63 

CTE- temperature exponent 

11 

FEX20 

F5.2 

65-69 

Reserved 

for 

future 

use 

Line 3: 








1 

EXPONENTS 

A9 

1-9 

Mnemonic 

name 


2 

FEX21 

F5.2 

11-15 

Reserved 

for 

future 

use 

3 

FEX22 

F5.2 

17-21 

Reserved 

for 

future 

use 

4 

FEX23 

F5.2 

23-27 

Reserved 

for 

future 

use 

5 

FEX24 

F5.2 

29-33 

Reserved 

for 

future 

use 

6 

FEX25 

F5.2 

35-39 

Reserved 

for 

future 

use 

7 

FEX26 

F5.2 

41-45 

Reserved 

for 

future 

use 

8 

FEX27 

F5.2 

47-51 

Reserved 

f or 

future 

use 

9 

FEX28 

F5.2 

53-57 

Reserved 

for 

future 

use 

10 

FEX29 

F5.2 

59-63 

Reserved 

for 

future 

use 

11 

FEX30 

F5.2 

65-69 

Reserved 

for 

future 

use 

Line 4: 








1 

EXPONENTS 

A9 

1-9 

Mnemonic 

name 


2 

FEX31 

F5.2 

11-15 

Reserved 

for 

future 

use 
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Remarks: 


1. The various exponents listed are used in the multifactor interaction 

equation (Section 1.6.1). For example FEX4 would be the exponent used to 
determine the temperature effect on the fiber modulus. 


t FEX4 


E 

T 

fll 

_Mf 

E 

T 

fllo 

Mf 


2. The values of the exponents listed above represent the default values. 


3. The locations of the currently used exponents are shown below for 
convenient reference . 




Exponent Field Name 

(i.e. FEX ) 

Property 

Stress Rate 

Stress 

Temperature 

Modulus 


2 

3 

4 

Strength 


5 

6 

7 

Poisson* s Ratio 


8 

9 

10 

Heat Conductivity 


14 

15 

16 

Coeff. of Thermal 

Exp. 

17 

18 

19 


All other locations are reserved for future use. 
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3 , 1*9 


Fiber Data Record: OVER 


Description : Defines the end of the fiber section in the databank 


Format: 


OVER END OF FIBER PROPERTIES 


Field No . Field Name Format Columns Description 
1 OVER. . . A28 1-28 Mnemonic name 

Remarks : 

1. This record is used to differentiate between the fiber and matrix sections 
of the databank. This record is only inserted after the last fiber entry. 
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3.2.1 


Matrix Data Record: 


CODE 


Description : Defines matrix code name 


Format and Example : 


+ 

CODE COMMENT 

+ 

TI15 Optional description of matrix 

+ 


Field No. 

Field Name 

Format 

Columns 

Description 

1 

CODE 

A4 

1-4 

Matrix code name 

2 

COMMENT 

A76 

5-80 

Optional description of matrix 


Remarks: 


1. The CODE for each matrix in the databank must be unique. 
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3.2.2 


Matrix Data Record: 


MP 


Description : Defines matrix physical properties 


Format and Example : 


-+ + 


MP p 

r m 

-+ + 

MP . 172E+00 
_+ + 


Field No. Field Name Format Columns 


1 FP A2 1-2 

3 D E10.3 4-13 


Description 

Mnemonic name 
Matrix weight density 
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3.2. 3 


Matrix Data Record: 


ME 


Description : Defines matrix elastic properties 


Format and Example : 


-+ 


ME 

-+ 


ME 

-+ 


+ — 

+ 

+ 

E 

V 

a 

m 

m 

m 




. 123E+08 
+ 

. 320E+00 
+ 

. 450E-05 

+ 


Field No. 

Field Name 

Format 

Columns 

Description 

1 

ME 

A2 

1-2 

Mnemonic name 

2 

E 

m 

E10.3 

4-13 

Matrix modulus 

3 

V 

m 

E10.3 

14-23 

Matrix Poisson's ratio 

4 

a 

m 

E10.3 

24-33 

Matrix coefficient of thermal 


expansion 






3.2.4 Matrix Data Record: 


Description: Defines matrix thermal properties 


Format and Example: 


m pm 

+ + 

. 390E+00 . 120E +00 


Field No. 

Field Name 

Format 

Columns 

Description 

1 

MT 

A2 

1-2 

Mnemonic name 

2 

K 

m 

E10.3 

4-13 

Matrix thermal conductivity 

3 

C 

pm 

E10.3 

14-23 

Matrix heat capacity 


MT 

-+ 

MT 
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3.2.5 Matrix Data Record: MS 


Description : Defines matrix strength properties 


Format and Example : 


1 

1 

+ 

1 H 

I i n B 

+ — 

S mC 

J U_ 

S mS 

+ 

e mT 

+ 

e mC 

+ 

G mS 

+ 

e 

mr 

; . 130E+06 

. 130E+06 

. 780E+05 

. 120E+00 

. 120E+00 

. 120E+00 

. 120E+00 

+ — 

+ 

+-- 

+ 

+ 

+ 

+ 


Field No. 

Field Name 

Format 

Columns 

Description 

1 

MS 

A2 

1-2 

Mnemonic name 

2 

S mT 

E10.3 

4-13 

Matrix 

tensile strength 

3 

S mC 

E10.3 

14-23 

Matrix 

compressive strength 

4 

S mS 

E10.3 

24-33 

Matrix 

shear strength 

5 

C mT 

E10.3 

34-43 

Matrix 

tensile strain 

6 

G mC 

E10.3 

44-53 

Matrix 

compressive strain 

7 

G mS 

E10.3 

54-63 

Matrix 

shear strain 

8 

e 

mx 

E10.3 

64-73 

Matrix 

torsional strain 


Remarks : 

1. In cases where only the experimental tensile strength is available, the 
following approximate values are used for the compressive and shear 
strengths: 

S c = S T and S g = 0.6 S T 
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3.2.6 


Matrix Data Record: 


MV 


Description : Defines additional matrix thermal properties 


Format and Example : 


-+- 


MV K 

v 


T 


Mm 


-+- 


+ 


MV . 019E+00 , 180E +04 

-+ + + 


Field No. 

Field Name 

Format 

Columns 

Description 

1 

MV 

A2 

1-2 

Mnemonic name 

2 

K 

V 

E10.3 

4-13 

Matrix void thermal conductivity 

3 

T 

Mm 

E10.3 

14-23 

Matrix melting temperature 
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3.2.7 Matrix Data Record: SIGMO 


Description : Defines reference values of matrix stresses 


Format and Example : 


SIGMO 

°mllA 

-+- 

°'m22A 

^228 

— + — 

°m2 2C 

-+ + — 

°ml2A 

+ — + - 

+ 

°'ml2B 

^120 

SIGMO 

0.0 


0.0 

0.0 


0.0 

0.0 

0.0 

0.0 

1 _L. _ -L L. 4. 

SIGMO 

‘’’n^A 


^238 

%23C 


^mlSA 

^mlSB 

_ .f. 

"mlSC 

°m33A 

SIGMO 

0.0 


0.0 

0.0 


0.0 

0.0 

4. — _ _ _ 4. . 

0.0 

0.0 


SIGMO 

‘’mOSB 


%i33C 

^ml IB 


°mllC 

+ — — 



SIGMO 

0.0 


0.0 

0.0 


0.0 

+ _ ^ . 













Field No. 

Field 

Name Format 

Columns 

Description 



Line 1 : 










1 

SIGMO 


A5 

1 

-5 

Mnemonic name 



2 

°'mllA 


E10. 3 

11 

-20 

Reference stress 

in region 

mllA 

3 

°'m22A 


E10. 3 

21 

-30 

Reference stress 

in region 

m22A 

4 

°m22B 


E10.3 

31 

-40 

Reference stress 

in region 

m22B 
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5 

%22C 

E10.3 

41-50 

Reference 

stress 

in 

region 

m22C 

6 

°ml2A 

E10.3 

51-60 

Reference 

stress 

in 

region 

ml2A 

7 

°ml2B 

E10.3 

61-70 

Ref erence 

stress 

in 

region 

ml2B 

8 

°'tnl2C 

E10.3 

71-80 

Reference 

stress 

in 

region 

ml2C 

Line 2: 









1 

SIGMO 

A5 

1-5 

Mnemonic name 




2 

°'m23A 

E10.3 

11-20 

Reference 

stress 

in 

region 

m23A 

3 

^238 

E10.3 

21-30 

Reference 

stress 

in 

region 

m23B 

4 

a 'm23C 

E10.3 

31-40 

Reference 

stress 

in 

region 

m23C 

5 

^mUA 

E10.3 

41-50 

Reference 

stress 

in 

region 

ml3A 

6 

‘VlBB 

E10.3 

51-60 

Reference 

stress 

in 

region 

ml3B 

7 

^mlBC 

E10.3 

61-70 

Reference 

stress 

in 

region 

ml3C 

8 

°m33A 

E10.3 

71-80 

Reference 

stress 

in 

region 

m33A 

Line 3: 









1 

SIGMO 

A5 

1-5 

Mnemonic name 




2 

^336 

E10.3 

11-20 

Reference 

stress 

in 

region 

m33B 

3 

%33C 

E10.3 

21-30 

Reference 

stress 

in 

region 

m33C 

4 

%11B 

E10.3 

31-40 

Reference 

stress 

in 

region 

ml IB 

5 

^mllC 

E10.3 

41-50 

Reference 

stress 

in 

region 

mllC 


Remarks: 


1. The reference values of stress are used in the stress term of the 

multifactor interaction relationship and are typically set to zero for 
convenience. 
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3.2,8 


Matrix Data Record: 


DOTMO 


Description : Defines reference values of matrix stress rates 


Format and Example : 


DOTMO 

VllA 

- + “ 

°m22A 

- + “ 

^mZZB 

- + - 

^10220 

- + - 

°ml2A 

*— + *— 

^mlZB 


%12C 















DOTMO 

0.0 


0.0 


0.0 


0.0 

X. 

0.0 


0.0 


0.0 



_ X — 








_ 


— * + — 


+ 


— T 


T 


T 








DOTMO 

^2 3 A 

. , 1. » 

^236 


^230 

.X*. 

^mlSA 

4. 

^mlSB 



^130 

4_ 

°m33A 

+ 




f 


T 


T 




T 


DOTMO 

0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


.+ + + + + + + 


DOTMO 

+ - 

^336 

^330 ^mllB 

- + ~ 

+ 

^mllC 

+ 

1 








DOTMO 

i 

O 1 

O I 

1 

0.0 0.0 
+ 

- + - 

0.0 

+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


Field No. Field Name Format 


Columns 


Description 


Line 1: 

1 DOTMO A5 


1-5 Mnemonic name 


^mllA 

^mZZA 


E10.3 11-20 Reference stress rate in region 

mllA 

E10.3 21-30 Reference stress rate in region 

m22A 
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°m22B 

°m22C 

^ml 2C 


E10.3 31-40 


E10.3 41-50 


E10.3 51-60 


E10.3 61-70 


E10.3 71-80 


Reference stress rate 
m22B 

Reference stress rate 
m22C 

Reference stress rate 
m!2A 

Reference stress rate 
ml2B 

Reference stress rate 
ml2C 


in region 


in region 


in region 


in region 


in region 


Line 2: 
1 
2 


3 


4 


5 


6 


7 


8 


DOTMO 

°'m23A 

^238 

°m23C 

^mlBA 

<r ml3B 

^mlBC 


A5 

E10.3 


1-5 

11-20 


E10.3 21-30 


E10.3 31-40 


E10.3 41-50 


E10.3 51-60 


E10.3 61-70 


Mnemonic name 

Reference stress rate in region 
m23A 

Reference stress rate in region 
m23B 

Reference stress rate in region 
m23C 

Reference stress rate in region 
ml3A 

Reference stress rate in region 
ml3B 

Reference stress rate in region 


ml3C 

E10.3 71-80 Reference stress rate in region 

m33A 
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Line 3: 


1 

DOTMO 

A5 

1-5 

Mnemonic name 




2 


E10. 3 

11-20 

Reference 

m33B 

stress 

rate 

in 

region 

3 

0 m33C 

E10. 3 

21-30 

Reference 

m33C 

stress 

rate 

in 

region 

4 

^mllB 

E10. 3 

31-40 

Reference 

mllB 

stress 

rate 

in region 

5 

cr mllC 

E10.3 

41-50 

Reference 
ml 1C 

stress 

rate 

in 

region 


Remarks: 


1. The reference values of stress rates are used in the stress rate term of 

the multifactor interaction relationship and are typically set to zero for 
convenience. 
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3 . 2.9 


Matrix Data Record: 


EXPONENTS 


Description : Defines values of exponents for matrix properties in the 

multifactor interaction relationship 

Format and Example : 


+- 

EXPONENTS 

MEX1 

+- 

MEX2 

+. 

MEX3 

MEX4 

+_ 

MEX5 

+- 

MEX6 

+ - 

MEX7 

MEX8 

+- 

MEX9 

+ 

MEX10 

EXPONENTS 

0.50 

0.50 

0.50 

0.50 

0.50 

0.00 

0.50 

0.50 

0.50 

0.50 

. 1 _ t j_ J J. 4- — 

EXPONENTS 
+ - 

MEX11 

MEX12 

MEX13 

MEX14 
+. 

MEX15 
+ - 

MEX16 

MEX17 
+ - 

MEX18 

MEX19 
+ - 

MEX20 


EXPONENTS 0.50 0.50 0.50 0.50 0.00 0.50 0.50 0.00 0.50 0.50 

+ + + + + + + + + + + 


+ - 

+- 

+_ 

+- 

+- 

+ - 

+. 

+. 

+- 

+. 

+ 

EXPONENTS 

MEX21 

MEX22 

MEX23 

MEX24 

MEX25 

MEX26 

MEX27 

MEX28 

MEX29 

MEX30 






.... i 



+ _ 

+ _ 


EXPONENTS 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

” T 

0.50 


+ + + + + + + + + + + 




EXPONENTS : 

MEX31 









h 


+ 

+ 









EXPONENTS 
+ - 

0.50 









+ 
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Field No. 

Field Name 

Format 

Columns 

Description 

Line 1 : 

1 

EXPONENTS 

A9 

1-9 

Mnemonic name 

2 

MEX1 

F5..2 

11-15 

Reserved for future use 

3 

MEX2 

F5.2 

17-21 

Modulus-stress rate exponent 

4 

MEX3 

F5.2 

23-27 

Modulus-stress exponent 

5 

MEX4 

F5.2 

29-33 

Modulus-temperature exponent 

6 

MEX5 

F5.2 

35-39 

Strength-stress rate exponent 

7 

MEX6 

F5.2 

41-45 

Strength-stress exponent 

8 

MEX7 

F5.2 

47-51 

Strength-temperature exponent 

9 

MEX8 

F5.2 

53-57 

Poisson’s ratio-stress rate 
exponent 

10 

MEX9 

F5.2 

59-63 

Poisson’s ratio-stress exponent 

11 

MEX10 

F5.2 

65-69 

Poisson’s ratio-temperature 
exponent 

Line 2: 

1 

EXPONENTS 

A9 

1-9 

Mnemonic name 

2 

MEX11 

F5.2 

11-15 

Reserved for future use 

3 

MEX12 

F5.2 

17-21 

Reserved for future use 

4 

MEX13 

F5.2 

23-27 

Reserved for future use 

5 

MEX14 

F5.2 

29-33 

Heat conductivity-stress rate 
exponent 

6 

MEX15 

F5.2 

35-39 

Heat conductivity-stress exponent 

7 

MEX16 

F5.2 

41-45 

Heat conduct ivity-temperature 
exponent 

8 

MEX17 

F5.2 

47-51 

CTE-stress rate exponent 
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9 

MEX18 

F5.2 

53-57 

CTE-stress rate exponent 

10 

MEX19 

F5.2 

59-63 

CTE-temperature exponent 

11 

MEX20 

F5.2 

65-69 

Reserved 

for future use 

Line 3: 

i 

EXPONENTS 

A9 

1-9 

Mnemonic 

name 

2 

MEX21 

F5.2 

11-15 

Reserved 

for future use 

3 

MEX22 

F5.2 

17-21 

Reserved 

for future use 

4 

MEX23 

F5.2 

23-27 

Reserved 

for future use 

5 

MEX24 

F5.2 

29-33 

Reserved 

for future use 

6 

MEX25 

F5.2 

35-39 

Reserved 

for future use 

7 

MEX26 

F5.2 

41-45 

Reserved 

for future use 

8 

MEX27 

F5.2 

47-51 

Reserved 

for future use 

9 

MEX28 

F5.2 

53-57 

Reserved 

for future use 

10 

MEX29 

F5.2 

59-63 

Reserved 

for future use 

11 

MEX30 

F5.2 

65-69 

Reserved 

for future use 

Line 4: 

1 

EXPONENTS 

A9 

1-9 

Mnemonic 

name 

2 

MEX31 

F5.2 

11-15 

Reserved 

for future use 
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Remarks: 


1. The various exponents listed are used in the multifactor interaction 

equation. For example MEX3 would be the exponent used to determine the 
stress effect on the matrix modulus. 

r n MEX3 


E 

s 

cr „ 

ml 1 

ml 1 

ml 1 

E 

s 

cr . . 

mil 

mil 

mil 


o o 


2. The values of the exponents listed above represent the default values. 

3. The locations of the currently used exponents are shown below for 
convenient ref erence . 




Exponent Field Name 

(i.e. MEX ) 

Property 

Stress Rate 

Stress 

Temperature 

Modulus 


2 

3 

4 

Strength 


5 

6 

7 

Poisson’s Ratio 


8 

9 

10 

Heat Conductivity 


14 

15 

16 

Coeff. of Thermal 

Exp. 

17 

18 

19 


All other locations are reserved for future use. 
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3 . 2.10 


Matrix Data Record: 


OVER 


Description : Defines the end of the matrix section in the databank 


Format: 


OVER END OF MATRIX PROPERTIES 


Field No . Field Name Format Columns Description 
1 OVER. . . A28 1-28 Mnemonic name 

Remarks : 

1. This record is used to differentiate between the matrix and interphase 
sections of the databank. This record is inserted only after the last 
matrix entry. 
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3.3.1 


Interphase Data Record: (No Mnemonic) 


Description : Defines optional interphase description 

Format and Example : 

+ + + + + + + + + + 

TITLE 

+ H + + — — — -4- H H — h 

Description of interphase 

+ + + + + + + + + + 

Field No. Field Name Format Columns Description 

1 TITLE A80 1-80 Interphase description 

Remarks : 

1. Unlike the fiber and matrix, which are distinguished by their code names, 
the interphases are defined by their thickness (see DP record). 
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3.3.2 Interphase Data Record: DP 


Description : Defines interphase physical properties 


Format and Example : 

-+ + + 

dp "a i d 

-+ + + 

DP . 172E+00 . 020E+00 

-+ + + 


Field No. 

Field Name 

Format 

Columns 

Description 

1 

DP 

A2 

1-2 

Mnemonic name 

2 

P d 

E10.3 

4-13 

Interphase weight density 

3 


E10.3 

14-23 

Interphase thickness 


Remarks : 

1. The input data file uses the interface thickness to select the desired 

interphase. This means that the thickness of the desired interphase must 
be unique. 
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3.3.3 Interphase Data Record: DE 


Description : Defines interphase elastic properties 


Format and Example : 


-+ — 




DE 

E d 

"d 

“d 

- + 




DE 
-+ — 

. 300E+07 

. 320E+00 
+ 

. 450E-05 


Field No. 

Field Name 

Format 

Columns 

Description 

1 

DE 

A2 

1-2 

Mnemonic name 

2 

E d 

E10.3 

4-13 

Interphase modulus 

3 

"d 

E10.3 

14-23 

Interphase poisson’s ratio 

4 

a d 

E10.3 

24-33 

Interphase coefficient of thermal 
expansion 
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Field No. 

Field Name 

Format 

Columns 

1 

DT 

A2 

1-2 

2 

K d 

E10.3 

4-13 

3 

C, 

E10.3 

14-23 


pd 
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Description 


Mnemonic name 

Interphase thermal conductivity 


Interphase heat capacity 




3.3.5 Interphase Data Record: DS 


Description : Defines interphase strength properties 


Format and Example : 


-+- 



- + + - 






+- 

+ 

DS 


S dT 

S dC 

S dS 


e dT 

L 

e dC 

C dS 

€ dx 

DS 


325E+05 . 325E+05 

. 195E+05 

. 120E+00 

. 120E+00 

. 120E+00 

. 120E+00 

+ 











Field 

No. 

Field Name 

Format 


Columns 


Description 



1 


DS 

A2 


1-2 


Mnemonic name 



2 


S dT 

E10.3 


4-13 


Interphase 

tensile 

strength 


3 


S dC 

E10.3 


14-23 


Interphase 

compressive strength 


4 


S dS 

E10.3 


24-33 


Interphase 

shear strength 


5 


e dT 

E10.3 


34-43 


Interphase 

tensile 

strain 


6 


e dC 

E10.3 


44-53 


Interphase 

compressive strain 


7 


e dS 

E10.3 


54-63 


Interphase 

shear strain 


8 


C dr 

E10.3 


64-73 


Interphase 

torsional strain 


Remarks : 

1. In cases where only the experimental tensile strength is available, the 
following approximate values are used for the compressive and shear 
strengths: 

S C * S T S s = 0 - 6S T 
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3.3.6 


Interphase Data Record: DV 


Description: Defines additional interphase thermal properties 


Format and Example: 


-+- 


+ 


DV K 

v 


T 


Md 


— +- 


+ 


DV . 019E+00 . 180E+04 

-+ + +■ 


Field No. 

Field Name 

Format 

Columns 

Description 

1 

DV 

A2 

1-2 

Mnemonic name 

2 

K 

V 

E10.3 

4-13 

Interphase void thermal 
conductivity 

3 

T 

Md 

E10.3 

14-23 

Interphase melting temperature 
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3 . 3.7 


Interphase Data Record: SIGDO 


Description : Defines reference values of interphase stresses 


Format and Example : 


SIGDO 

<r dllB °d22B 

^220 

- + - 

°‘dl2B 

- + + + — 

°dl2C °d23B 

^230 

+ 




T 






SIGDO 

0.0 0. 

,0 

0.0 

-+- 

0.0 

0. 

“ + 

O 1 

1 
1 
1 

+ 

1 

1 

1 

O 1 
1 

,0 

0.0 


+ + + + + + + 


SIGDO 

°dl3B 

^130 

+ + 

°d33B 

^330 "dllC 













SIGDO 

0.0 

0.0 

0.0 

0.0 0. 

,0 













Field No. 

Field 

Name Format 

Columns 

Description 




Line Is 









1 

SIGDO 

A5 

1-5 

Mnemonic name 




2 

°dllB 

E10.3 

11-20 

Reference 

stress 

in 

region 

dllB 

3 

°d22B 

E10.3 

21-30 

Reference 

stress 

in 

region 

d22B 

4 

°d22C 

ElU. 3 

31-40 

Reference 

stress 

in 

region 

d22C 

5 

‘’’d^B 

E10.3 

41-50 

Reference 

stress 

in 

region 

dl2B 

6 

^120 

E10.3 

51-60 

Reference 

stress 

in 

region 

dl2C 

7 

°d23B 

E10.3 

61-70 

Reference 

stress 

in 

region 

d23B 

8 

^230 

E10.3 

71-80 

Reference 

stress 

in 

region 

d23C 


+ 


+ 


+ 


+ 


+ 
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Line 2: 


1 

SIGDO 

A5 

1-5 

Mnemonic name 



2 

^dlSB 

E10.3 

11-20 

Reference 

stress 

in region 

dl3B 

3 

°dl3C 

E10.3 

21-30 

Reference 

stress 

in region dl3C 

4 

°d33B 

E10.3 

31-40 

Reference 

stress 

in region 

d33B 

5 

°d33C 

E10.3 

41-50 

Reference 

stress 

in region 

d33C 

6 

^dllC 

E10.3 

51-60 

Reference 

stress 

in region 

dl 1C 


Remarks : 

^ Tiig reference values of the stresses are used in the stress term of the 
multifactor interaction relationship and are typically set to zero for 
convenience. 
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3.3.8 Interphase Data Record: DOTDO 


Description : Defines reference values of interphase stress rates 


Format and Example : 


DOTDO 

°dl IB 

- + - 
— -f — 

+ _ 

^226 

°"d22C 

- + - 

0* 

dl2B 

- + 

& 

dl2C 

- + - 

°d23B 

-+- 

°d23C 

DOTDO 

0.0 

- + - 

0.0 

+_ 

0.0 

r 

- + - 

0.0 

0.0 

- + 

— + — 
- + - 

0.0 

— + — 
— h— 

0.0 


+ + + + + + 


DOTDO 

^dlBB 

°'dl3C 

°d33B 

°d33C "dllC 



DOTDO 

0.0 

0.0 

— -I + _ 

0.0 

0.0 0.0 



Field No. 

Field 

Name Format 

Columns 

Description 



Line 1 : 







1 

DOTDO 

A5 

1-5 

Mnemonic name 



2 

^dllB 

E10.3 

11-20 

Reference stress 

rate 

in region 





dllB 



3 

°d22B 

E10.3 

21-30 

Reference stress 

rate 

in region 





d22B 



4 

°d22C 

E10.3 

31-40 

Reference stress 

rate 

in region 





d22C 



5 

°dl2B 

E10.3 

41-50 

Reference stress 

rate 

in region 





d!2B 
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6 *dl2C E10 ’ 3 51-60 

7 *d23B E1 °- 3 61 " 7 ° 

8 *d23C E1 °- 3 71-80 


Reference stress rate in region 
dl2C 

Reference stress rate in region 
d23B 

Reference stress rate in region 
d23C 



A5 1-5 Mnemonic name 

E10.3 11-20 Reference stress rate in region 

dl3B 

E10.3 21-30 Reference stress rate in region 

dl3C 

E10.3 31-40 Reference stress rate in region 

d33B 

E10.3 41-50 Reference stress rate in region 

d33C 

E10.3 51-60 Reference stress rate in region 

dllC 


1. The reference values of the stress rates are used in the stress rate term 
of the multifactor interaction relationship and are typically set to zero 
for convenience. 
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3.3.9 Interphase Data Record: EXPONENTS 


Description : Defines values of exponents for interphase properties in 

multifactor interaction relationship 


Format and Example: 


EXPONENTS DEX1 DEX2 DEX3 DEX4 DEX5 DEX6 DEX7 DEX8 DEX9 DEX10 

+ + + + + + + + + + + 

EXPONENTS 0.50 0.50 0.50 0.50 0.50 0.00 0.50 0.50 0.50 0.50 



EXPONENTS DEX11 DEX12 DEX13 DEX14 DEX15 DEX16 DEX17 DEX18 DEX19 DEX20 
+ + + + + + + + + + + 

EXPONENTS 0.50 0.50 0.50 0.50 0.00 0.50 0.50 0.00 0.50 0.50 



EXPONENTS 

DEX21 

DEX22 

1- 

DEX23 

- - - -L 

DEX24 

_L 

DEX25 

DEX26 

DEX27 

DEX28 

DEX29 

DEX30 

EXPONENTS 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

■ “ “ “ ™ X* 

0.50 

0.50 

0.50 
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Field No. 

Field Name 

Format 

Columns 

Description 

Line 1: 

1 

EXPONENTS 

A9 

1-9 

Mnemonic name 

2 

DEX1 

F5.2 

11-15 

Reserved for future use 

3 

DEX2 

F5.2 

17-21 

Modulus-stress rate exponent 

4 

DEX3 

F5.2 

23-27 

Modulus-stress exponent 

5 

DEX4 

F5.2 

29-33 

Modulus-temperature exponent 

6 

DEX5 

F5.2 

35-39 

Strength-stress rate exponent 

7 

DEX6 

F5.2 

41-45 

Strength-stress exponent 

8 

DEX7 

F5.2 

47-51 

Strength-temperature exponent 

9 

DEX8 

F5.2 

53-57 

Poisson's ratio-stress rate 
exponent 

10 

DEX9 

F5.2 

59-63 

Poisson's ratio-stress exponent 

11 

DEX10 

F5.2 

65-69 

Poisson's ratio-temperature 
exponent 

Line 2: 

1 

EXPONENTS 

A9 

1-9 

Mnemonic name 

2 

DEX11 

F5.2 

11-15 

Reserved for future use 

3 

DEX12 

F5.2 

17-21 

Reserved for future use 

4 

DEX13 

F5.2 

23-27 

Reserved for future use 

5 

DEX14 

F5.2 

29-33 

Heat conductivity-stress rate 
exponent 

6 

DEX15 

F5.2 

35-39 

Heat conductivity-stress exponent 

7 

DEX16 

F5.2 

41-45 

Heat conductivity-temperature 
exponent 

8 

DEX17 

F5.2 

47-51 

CTE-stress rate exponent 
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9 

DEX18 

F5.2 

53-57 

CTE-stress rate exponent 

10 

DEX19 

F5.2 

59-63 

CTE-temperature exponent 

11 

DEX20 

F5.2 

65-69 

Reserved 

for future 

use 

Line 3: 







1 

EXPONENTS 

A9 

1-9 

Mnemonic 

name 


2 

DEX21 

F5.2 

11-15 

Reserved 

for future 

use 

3 

DEX22 

F5.2 

17-21 

Reserved 

for future 

use 

4 

DEX23 

F5.2 

23-27 

Reserved 

for future 

use 

5 

DEX24 

F5.2 

29-33 

Reserved 

for future 

use 

6 

DEX25 

F5.2 

35-39 

Reserved 

for future 

use 

7 

DEX26 

F5.2 

41-45 

Reserved 

for future 

use 

8 

DEX27 

F5.2 

47-51 

Reserved 

for future 

use 

9 

DEX28 

F5.2 

53-57 

Reserved 

for future 

use 

10 

DEX29 

F5.2 

59-63 

Reserved 

for future 

use 

11 

DEX30 

F5.2 

65-69 

Reserved 

for future 

use 

Line 4: 







1 

EXPONENTS 

A9 

1-9 

Mnemonic 

name 


2 

DEX31 

F5.2 

11-15 

Reserved 

for future 

use 
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Remarks: 


1. The various exponents listed are used in the multifactor interaction 

equation. For example DEX4 would be the exponent used to determine the 
stress effect on the matrix modulus. 


"dll 


"dll 


n DEX3 


S dll" <r dll 


L S dll“ °dll -> 
o 


2. The values of the exponents listed above represent the default values. 


3. The locations of the currently used exponents are shown below for 
convenient reference. 




Exponent 

Field Name 

(i.e. DEX ) 

Property 

Stress Rate 

Stress 

Temperature 

Modulus 


2 

3 

4 

Strength 


5 

6 

7 

Poisson’s Ratio 


8 

9 

10 

Heat Conductivity 


14 

15 

16 

Coeff. of Thermal 

Exp. 

17 

18 

19 


All other locations are reserved for future use. 


3-47 



3.4 Sample Constituent Databank 


A sample databank is provided in Figure 3.4-1 with entries for two fibers (SiC 
and Tungsten), two matrices (Ti-15-3 and Copper), and an interphase (a carbon 
coating). The databank begins with the fiber properties, followed by the matrix 
properties, and ends with the interphase properties. 
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SICA SILICON CARBIDE FIBER 


FP 1 0 . 560E-02 O.llOE+OO 0.487E+04 

FE 0. 620E+08 0.620E+08 0.300E+00 0.300E+00 0.238E+08 0.238E+08 

FT 0. 272E-05 0.272E-0S 0.750E+00 0.750E+00 0.290E+00 

FS 0. 500E+06 0 . 650E+06 0.500E+06 0.650E+06 0.300E+06 0.300E+06 


SIGFO 


0.0 

0.0 

0.0 


0.0 

0.0 

0.0 

DOTFO 


0.0 

0.0 

0.0 


0.0 

0.0 

0.0 

EXPONENTS 

0.25 

0.25 

0.25 

0.25 

0.25 

0.00 

0.25 

0.25 

0.25 

0.25 

EXPONENTS 

0.25 

0.25 

0.25 

0.25 

0.00 

0.25 

0.25 

0.00 

0.25 

0.25 

EXPONENTS 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 


EXPONENTS 0.25 
TUNG TUNGSTEN FIBER 

FP 1 0. 100E-01 0. 683E+00 0.617E+04 

FE 0. 590E+08 0.590E+08 0.290E+00 0.290E+00 0.227E+08 0.227E+08 

FT 0. 250E-05 0.250E-05 0.828E+01 0.828E+01 0.240E-01 

FS 0. 390E+06 0.390E+06 0.390E+06 0.390E+06 0.236E+06 0.236E+06 


SIGFO 


0.0 

0.0 

0.0 


0.0 

0.0 

0.0 

DOTFO 


0.0 

0.0 

0.0 


0.0 

0.0 

0.0 

EXPONENTS 

0.25 

0.25 

0.25 0.25 

0.25 

0.00 

0.25 

0.25 0.25 

0.25 

EXPONENTS 

0.25 

0.25 

0.25 0.25 

0.00 

0.25 

0.25 

0.00 0.25 

0.25 

EXPONENTS 

0.25 

0.25 

0.25 0.25 

0.25 

0.25 

0.25 

0.25 0.25 

0.25 

EXPONENTS 

0.25 









OVER END OF FIBER PROPERTIES 
TI15 TITANIUM MATRIX 
MP 0. 170E+00 

ME 0 . 1 23E+Q8 0.320E+00 0.450E-05 
MT 0 . 390E+00 0.120E+00 


MS 0 . 130E+06 

0. 130E+06 0. 

. 780E+05 0. 

120E+00 0. 

120E+00 0. 

120E+00 0. 

120E+00 


MV 0 . 019E+00 

0. 180E+04 







SIGMO 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

SIGMO 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

SIGMO 

0.0 

0.0 

0.0 

0.0 




DOTMO 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

DOTMO 

0.0 

0.0 

0.0 

0.0 

0.0 

o 

o 

0.0 

DOTMO 

0.0 

0.0 

0.0 

0.0 





Figure 3.4-1 — Sample case constituent databank. 
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EXPONENTS 0.50 0.50 0.50 0.50 0.50 0.00 0.50 0.50 0.50 0.50 

EXPONENTS 0.50 0.50 0.50 0.50 0.00 0.50 0.50 0.00 0.50 0.50 

EXPONENTS 0.50 0.50 0.50 0.50 0.50 0.00 0.50 0.50 0.50 0.50 


EXPONENTS 0.50 
COPR COPPER MATRIX 
HP 0. 320E+00 

ME 0. 177E+08 0.300E+00 0.980E-05 
MT 0. 193E+02 0.090E+00 

MS 0. 320E+05 0.320E+05 0.190E+05 0.350E+00 0.350E+00 0.350E+00 0.350E+00 
MV 0. 019E+00 0. 198E+04 


SIGMO 


0.0 

0.0 

0.0 


0.0 

0.0 

0.0 

0.0 

SIGMO 


0.0 

0.0 

0.0 


0.0 

0.0 

0.0 

0.0 

SIGMO 


0.0 

0.0 

0.0 


0.0 




DOTMO 


0.0 

0.0 

0.0 


0.0 

0.0 

0.0 

0.0 

DOTMO 


0.0 

0.0 

0.0 


0.0 

0.0 

0.0 

0.0 

DOTMO 


0.0 

0.0 

0.0 


0.0 




EXPONENTS 

0.50 

0.50 

0.50 

0.50 

0.50 

0.00 

0.50 

0.50 0.50 

0. 50 


EXPONENTS 

0.50 

0.50 

0.50 

0.50 

0.00 

0.50 

0.50 

0.00 0.50 

0.50 


EXPONENTS 

0.50 

0.50 

0.50 

0.50 

0.50 

0.00 

0.50 

0.50 0.50 

0.50 



EXPONENTS 0.50 

OVER END OF MATRIX PROPERTIES 

INTERPHASE CARBON COATING FOR FIBERS 

DP 0 . 172E+00 0 . 020E+0Q 

DE 0 . 250E+07 0.220E+00 0.212E-05 

DT 0. 390E+00 0.120E+00 

DS 0. 100E+05 0. 100E+05 0. 100E+05 0. 120E+00 0. 120E+00 0. 120E+00 0. 120E+00 
DV 0.019E+00 0. 140E+04 


SIGDO 


0.0 

0.0 

0.0 


0.0 

0.0 

0.0 

SIGDO 


0.0 

0.0 

0.0 


0.0 

0.0 

0.0 

DQTDO 


0.0 

0.0 

0.0 


0.0 

0.0 

0.0 

DOTDO 


0.0 

0.0 

0.0 


0.0 

0.0 

0.0 

EXPONENTS 

0.50 

0.50 

0.50 0.50 

0.50 

0.00 

0.50 

0.50 0.50 

0.50 

EXPONENTS 

0.50 

0.50 

0.50 0.50 

0.00 

0.50 

0.50 

0.00 0.50 

0.50 

EXPONENTS 

0.50 

0.50 

0.50 0.50 

0.50 

0.00 

0.50 

0.50 0.50 

0.50 

EXPONENTS 

0.50 









Figure 3.4-1 Concluded. 
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4.0 Postprocessing Output Files 


The ten postprocessing files contain information pertaining to various 
behavior/response variables. As noted earlier, these files can be optionally 
generated at the user’s convenience and essentially duplicates the same 
information which is reported in the primary output file, but in format from 
which selected output quantities can be extracted by the user and imported into 
graphics packages for convenient representation. The general format of each 
line of the postprocessing files is ( 5 ( IX , El 5 . 6 ) ) . The first entry in each of 
the ten files contains the time for the current step. A description of the 
contents of each file is given below. 
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4.1 Postprocessing File 1 


Description : Time and temperature for each ply 


Format : 



fc i 

T 

ply i 

T 

ply 2 

T . 

piy n 

Field 

No. 

Field Name Format 

Columns 

Description 

1 


t i E15.6 

2-16 

Time at each load/time step 

2 


T . . E15.6 

ply 1 

18-32 

Temperature of ply 1 at time 

3 


T . _ E15.6 

ply 2 

34-48 

Temperature of ply 2 at time 

n 


T , E15.6 

ply n 


Temperature of ply n at time 


Remarks : 

1. If there are more than four plies, the temperatures for plies five onward 
are continued on subsequent lines until all plies are represented. 

2. The set of lines is repeated for each time step in the analysis. 
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4.2 Postprocessing File 2 


Description: Time and instantaneous fiber properties for each ply 


Format : 



Line 1: 


Field No. 

Field Name 

Format 

Columns 

1 

t. 

i 

E15.6 

2-16 

2 

27 

E15.6 

18-32 

3 

P f 

E15.6 

34-48 

4 

E fll 

E15.6 

50-64 

5 

E f22 

E15.6 

66-80 

Line 2: 




Field No. 

Field Name 

Format 

Columns 

1 

G f 12 

E15.6 

2-16 

2 

G f 23 

E15.6 

18-32 

3 

V f 12 

E15.6 

34-48 

4 

V f23 

E15.6 

50-64 

5 

C ff 

E15.6 

66-80 

Line 3: 




Field No. 

Field Name 

Format 

Columns 

1 

K fn 

E15.6 

2-16 

2 

K f22 

E15.6 

18-32 

3 

“fll 

E15.6 

34-48 

4 

“f22 

E15. 6 

50-64 

5 

S fllT 

E15.6 

66-80 


Description 

Time at each load/time step 
Number of fiber properties 
Fiber density at time t^ 

Fiber 11 modulus at time t. 

i 

Fiber 22 modulus at time t^ 
Description 

Fiber 12 shear modulus at time t. 

i 

Fiber 23 shear modulus at time t. 

i 

Fiber 12 Poisson's ratio at time t. 

l 

Fiber 23 Poisson's ratio at time t. 

l 

Fiber heat capacity at time t^ 
Description 

Fiber 11 thermal conductivity at 
time t. 

l 

Fiber 22 thermal conductivity at 
time t. 

i 

Fiber 11 coefficient of thermal 
expansion at time t^ 

Fiber 22 coefficient of thermal 
expansion at time t^ 

Fiber 11 tensile strength at 
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Line 4: 


Field No. Field Name Format Columns 


1 

S fllC 

E15.6 

2-16 

2 

S f22T 

E15.6 

18-32 

3 

S f22C 

E15.6 

34-48 

4 

S fl2 

E15.6 

50-64 

5 

S f23 

E15.6 

66-80 

Line 5: 




Field No. 

Field Name 

Format 

Columns 

1 

S fl3 

E15.6 

2-16 

2 

13 

E15.6 

18-32 

3 

G f 13 

E15.6 

34-48 

4 

E 

1 33 

E15.6 

50-64 

5 

a f 33 

E15.6 

66-80 


Line 6: 




Field No. 

Field Name 

Format 

Columns 

1 

G f 12 

E15.6 

2-16 

2 

G f23 

E15.6 

18-32 

3 

l fl3 

E15.6 

34-48 

4 

S f33T 

E15.6 

50-64 


4-6 


Fiber 33 tensile strength at time 



5 


E15.6 


66-80 


S f 33C 


Fiber 33 compressive strength at 
time t. 

l 


Remarks : 

1. The set of six lines is repeated for each time step in the analysis. 
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4.3 Postprocessing File 3 


Description: Time and instantaneous matrix properties for each ply 


Format: 





4-8 


G ml3C G ml3 V ml3 V ml2A "m^B 



V ml2C V m23A V m23B U m23C %13A 



E mllA E mllB E mllC S m33T S m33C 


Line 1: 


Field No. 

Field Name 

Format 

Columns 

Description 

1 

t. 

i 

E15.6 

2-16 

Time at each load/time step 

2 

59 

E15.6 

18-32 

Number of matrix properties 

3 

p m 

E15.6 

34-48 

Matrix density at time t. 

4 

E mll 

E15.6 

50-64 

Matrix 11 modulus at time t ^ 

5 

E m22A 

E15.6 

0^ 

ON 

1 

00 

o 

Matrix 22A modulus at time t. 

l 
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Line 2: 


Field No. 

Field Name 

Format 

Columns 

Description 

1 

E m22B 

E15.6 

2-16 

Matrix 22B modulus at time t^ 

2 

E m22C 

E15.6 

18-32 

Matrix 22C modulus at time t^ 

3 

E m22 

E15.6 

34-48 

Matrix 22 modulus at time t 

4 

G ml2A 

E15.6 

50-64 

Matrix 12A shear modulus at time t 

l 

5 

G ml2B 

E15.6 

66-80 

Matrix 12B shear modulus at time t 

i 

Line 3: 

Field No. 

Field Name 

Format 

Columns 

Description 

1 

G ml2C 

E15.6 

2-16 

Matrix 12C shear modulus at time t^ 

2 

G ml2 

E15.6 

18-32 

Matrix 12 shear modulus at time t. 

1 

3 

G m23A 

E15.6 

34-48 

Matrix 23A shear modulus at time 

4 

G m23B 

E15.6 

50-64 

Matrix 23B shear modulus at time 

5 

G m23C 

E15.6 

66-80 

Matrix 23C shear modulus at time t^ 

Line 4: 

Field No. 

Field Name 

Format 

Columns 

Description 

1 

G m23 

E15.6 

2-16 

Matrix 23 shear modulus at time t 

i 

2 

l ml2 

E15.6 

18-32 

Matrix 12 Poisson’s ratio at 
time t. 

i 

3 

V m23 

E15.6 

34-48 

Matrix 23 Poisson’s ratio at 
time t^ 

4 

C fm 

E15.6 

50-64 

Matrix heat capacity at time t^ 

5 

K mll 

E15.6 

66-80 

Matrix 11 thermal conductivity at 
time t^ 
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Line 5: 


Field No. 

Field Name 

Format 

Columns 

Description 

1 

K m22 

E15.6 

2-16 

Matrix 22 thermal conductivity at 
time t^ 

2 

“mil 

E15.6 

18-32 

Matrix 11 coefficient of thermal 
expansion at time t^ 

3 

“m22A 

E15.6 

34-48 

Matrix 22A coefficient of thermal 
expansion at time t^ 

4 

%22B 

E15.6 

50-64 

Matrix 22B coefficient of thermal 
expansion at time t ^ 

5 

a m22C 

E15.6 

66-80 

Matrix 22C coefficient of thermal 
expansion at time t^ 

Line 6: 

Field No. 

Field Name 

Format 

Columns 

Description 

1 

“m22 

E15.6 

2-16 

Matrix 22 coefficient of thermal 
expansion at time t^ 

2 

S ml IT 

E15.6 

18-32 

Matrix 11 tensile strength at 

time t. 

l 

3 

S mllC 

E15.6 

34-48 

Matrix 11 compressive strength at 
time t. 

i 

A 

-t. 

S m22T 

E15.6 

50-64 

Matrix 22 tensile strength at 
time t. 

l 

5 

S m22C 

E15.6 

66-80 

Matrix 22 compressive strength at 
time t^ 
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Line 7: 


Field No. 

Field Name 

Format 

Columns 

1 

S ml2 

E15.6 

2-16 

2 

S m23 

E15.6 

18-32 

3 

S ml3 

E15.6 

34-48 

4 

G ml3A 

E15.6 

50-64 

5 

G ml3B 

E15.6 

66-80 

Line 8: 

Field No. 

Field Name 

Format 

Columns 

1 

G ml3C 

E15.6 

2-16 

2 

G ml3 

E15.6 

18-32 

3 

t ’ml3 

E15.6 

34-48 

4 

l ml2A 

E15.6 

50-64 

5 

l 'ml2B 

E15.6 

66-80 


Description 


Matrix 

12 

shear 

strength 

at 

time 


Matrix 

23 

shear 

strength 

at 

time 

t. 

i 

Matrix 

13 

shear 

strength 

at 

time 

t. 

i 

Matrix 

13A 

shear 

modulus 

at 

time 

fc i 

Matrix 

13B 

shear 

modulus 

at 

time 

t. 


i 


Description 

Matrix 13C shear modulus at time t 
Matrix 13 shear modulus at time t^ 
Matrix 13 Poisson's ratio at 
time t^ 

Matrix 12A Poisson’s ratio at 
time t^ 

Matrix 12B Poisson’s ratio at 
time t. 

i 
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Line 9: 


Field No. 

Field Name 

Format 

Columns 

1 

V mlZC 

E15.6 

2-16 

2 

V m23A 

E15.6 

18-32 

3 

V m23B 

E15.6 

34-48 

4 

%23C 

E15.6 

50-64 

5 

y ml3A 

E15.6 

66-80 


Description 

Matrix 12C Poisson's ratio at 
time t^ 

Matrix 23A Poisson's ratio at 

time t. 

1 

Matrix 23B Poisson's ratio at 
time 

Matrix 23C Poisson's ratio at 

time t. 

i 

Matrix 13A Poisson's ratio at 
time t. 

i 


Line 10: 


Field No. 

Field Name 

Format 

Columns 

Description 



1 

‘'mlSB 

E15.6 

2-16 

Matrix 13B 

Poisson' 

s ratio 

at 





time t^ 




2 

^130 

E15.6 

18-32 

Matrix 13C 

Poisson* 

s ratio 

at 





time t. 
i 




3 

E m33A 

E15.6 

34-48 

Matrix 33A 

modulus 

at time 

t. 

i 

4 

E m33B 

E15.6 

50-64 

Matrix 33B 

modulus 

at time 

t . 

l 

5 

E m33C 

E15.6 

66-80 

Matrix 33C 

modulus 

at time 

t. 

l 
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Line 11: 


Field No 

Field Name 

Format 

Columns 

Description 

i 

E m33 

E15.6 

2-16 

Matrix 33 modulus at time t^ 

2 

a m33A 

E15.6 

18-32 

Matrix 33A coefficient of thermal 
expansion at time 

3 

a m33B 

E15.6 

34-48 

Matrix 33B coefficient of thermal 
expansion at time t^ 

4 

“m33C 

E15.6 

50-64 

Matrix 33C coefficient of thermal 
expansion at time t^ 

5 

“m33 

E15.6 

66-80 

Matrix 33 coefficient of thermal 
expansion at time t^ 

Line 12: 

Field No. 

Field Name 

Format 

Columns 

Description 

1 

E mllA 

E15.6 

2-16 

Matrix 11A modulus at time t^ 

2 

E mllB 

E15.6 

18-32 

Matrix 11B modulus at time t^ 

3 

E mllC 

E15.6 

34-48 

Matrix 11C modulus at time t^ 

4 

S m33T 

E15.6 

50-64 

Matrix 33 tensile strength at 
time t^ 

5 

S m33C 

E15.6 

66-80 

Matrix 33 compressive strength at 
time t^ 

Remarks : 

1. The 

set of twelve 

lines is 

repeated for 

each time step in the analysis. 
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4.4 Postprocessing File 4 


Description : Time and instantaneous interphase properties for each ply 


Format : 



4-15 









+ 


+ 


+ 


+ 


+ 


G G v v v 

dl3C dl3 dl3 dl2B d!2C 

+ + + + + 

+ + + + + 


^dZSB y dZ3C ^dlZB ‘'dISC E d33B 

4 + + + + 

+ + + + + 


E d33C E d33 “d33B “d33C “d33 

4 + + + + 

+ 4 4. 4 + 


E dllB 


E dllC 


S d33T 


S d33C 


+ 


+ 


Line 1: 
Field No. 

Field Name 

Format 

Columns 

Description 

1 

t. 

1 

E15.6 

Z-16 

Time at each load/time step 

Z 

48 

E15.6 

18-3Z 

Number of interphase properties 

3 

p d 

E15.6 

34-48 

Interphase density at time t^ 

4 

E dll 

E15.6 

50-64 

Interphase 11 modulus at time t^ 

5 

E d2ZB 

E15.6 

66-80 

Interphase 22B modulus at time t^ 
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Line 2: 


Field No. 

Field Name 

Format 

Columns 

Description 

1 

E d22C 

E15.6 

2-16 

Interphase 22C modulus at time t^ 

2 

E d22 

E15.6 

18-32 

Interphase 22 modulus at time t^ 

3 

G dl2B 

E15.6 

34-48 

Interphase 12B shear modulus at 

time t. 

l 

4 

G dl2C 

E15.6 

50-64 

Interphase 12C shear modulus at 
time 

5 

Line 3: 

G dl2 

E15.6 

66-80 

Interphase 12 shear modulus at 
time t. 

l 

Field No. 

Field Name 

Format 

Columns 

Description 

1 

G d23B 

E15.6 

2-16 

Interphase 23B shear modulus at 
time t. 

l 

2 

G d23C 

E15.6 

18-32 

Interphase 23C shear modulus at 
time tj. 

3 

G d23 

E15.6 

34-48 

Interphase 23 shear modulus at 
time t. 

l 

4 

P dl2 

E15.6 

50-64 

Interphase 12 Poisson's ratio at 
time t^ 

5 

y d23 

E15.6 

66-80 

Interphase 23 Poisson’s ratio at 
time t. 

i 
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Line 4: 


Field No. 

Field Name 

Format 

Columns 

1 

c fd 

E15.6 

2-16 

2 

K dll 

E15.6 

18-32 

3 

K d22 

E15.6 

34-48 

4 

“dll 

E15.6 

50-64 

5 

“d22B 

E15.6 

66-80 

Line 5: 




Field No. 

Field Name 

Format 

Columns 

i 

“d22C 

E15.6 

2-16 

2 

“d22 

E15.6 

18-32 

3 

S dllT 

E15.6 

34-48 

4 

S dl 1C 

E15.6 

50-64 

5 

S d22T 

E15.6 

66-80 


Description 

Interphase heat capacity at time t^ 
Interphase 11 thermal conductivity 
at time t^ 

Interphase 22 thermal conductivity 
at time t^ 

Interphase 11 coefficient of 
thermal expansion at time t^ 
Interphase 22B coefficient of 
thermal expansion at time t^ 

Description 

Interphase 22C coefficient of 
thermal expansion at time t^ 
Interphase 22 coefficient of 
thermal expansion at time t^ 
Interphase 11 tensile strength at 
time t^ 

Interphase 11 compressive strength 
at time t^ 

Interphase 22 tensile strength at 
time t^ 
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Line 6: 


Field No. 

Field Name 

Format 

Columns 

Description 

1 

S d22C 

E15.6 

2-16 

Interphase 
at time t^ 

22 compressive strength 

2 

S dl2 

E15.6 

18-32 

Interphase 
time t^ 

12 shear strength at 

3 

S d23 

E15.6 

34-48 

Interphase 
time t^ 

23 shear strength at 

4 

S dl3 

E15.6 

50-64 

Interphase 
time t^ 

13 shear strength at 

5 

Line 7: 

G dl3B 

E15.6 

66-80 

Interphase 
time t. 

i 

13B shear modulus at 

Field No. 

Field Name 

Format 

Columns 

Description 

1 

G dl3C 

E15.6 

2-16 

Interphase 
time t. 

l 

13C shear modulus at 

2 

G dl3 

E15.6 

18-32 

Interphase 
time t^ 

13 shear modulus at 

3 

V dl3 

E15.6 

34-48 

Interphase 
time t^ 

13 Poisson* s ratio at 

4 

^dl2B 

E15.6 

50-64 

Interphase 
time t. 

l 

12B Poisson* s ratio at 

5 

V 

d!2C 

E15.6 

66-80 

Interphase 
time t^ 

12C Poisson’ s ratio at 
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Line 8: 


Field No. 

Field Name 

Format 

Columns 

i 

V d23B 

E15.6 

2-16 

2 

V d23C 

E15.6 

18-32 

3 

^dl3B 

E15.6 

34-48 

4 

^dl3C 

E15.6 

50-64 

5 

E d33B 

E15.6 

66-80 

Line 9: 

Field No. 

Field Name 

Format 

Columns 

1 

E d33C 

E15.6 

2-16 

2 

E 

d33 

E15.6 

18-32 

3 

“d33B 

E15.6 

34-48 

4 

“d33C 

E15.6 

50-64 

5 

“d33 

E15.6 

66-80 


Description 

Interphase 23B Poisson’s ratio at 
time t^ 

Interphase 23C Poisson’s ratio at 
time t ^ 

Interphase 13B Poisson’s ratio at 
time t^ 

Interphase 13C Poisson’s ratio at 
time t^ 

Interphase 33B modulus at time t 


Description 

Interphase 33C modulus at time t^ 
Interphase 33 modulus at time t ^ 
Interphase 33B coefficient of 
thermal expansion at time t^ 
Interphase 33C coefficient of 
thermal expansion at time t^ 
Interphase 33 coefficient of 
thermal expansion at time t^ 
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Line 10: 


Field No. 

Field Name 

Format 

Columns 

1 

E dllB 

E15.6 

2-16 

2 

E dllC 

E15.6 

18-32 

3 

S d33T 

E15.6 

34-48 

4 

S d33C 

E15.6 

50-64 

Remarks: 




1. The 

set of ten lines is repeated for 


Description 

Interphase 11B modulus at time t^ 
Interphase 11C modulus at time t ^ 
Interphase 33 tensile strength at 
time t^ 

Interphase 33 compressive strength 
at time t^ 

each time step in the analysis. 
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4.5 Postprocessing File 5 


Description : Time and instantaneous ply properties for each ply 


Format: 



4-22 





Line 1 : 


Field No. 

Field Name 

Format 

Columns 

Description 

1 

fc i 

E15.6 

2-16 

Time at each load/time step 

2 

27 

E15.6 

18-32 

Number of ply properties 

3 

p l 

E15.6 

34-48 

Ply density at time t^ 

4 

E m 

E15.6 

50-64 

Ply 11 modulus at time t^ 

5 

E 122 

E15.6 

66-80 

Ply 22 modulus at time t ^ 

Line 2: 

Field No. 

Field Name 

Format 

Columns 

Description 

1 

E 133 

E15.6 

2-16 

Ply 33 modulus at time t^ 

2 

G 112 

E15.6 

18-32 

Ply 12 shear modulus at time t ^ 

3 

G 123 

E15.6 

34-48 

Ply 23 shear modulus at time t ^ 

4 

G 113 

E15.6 

50-64 

Ply 13 shear modulus at time t 

5 

V U2 

E15.6 

66-80 

Ply 12 Poisson' s* rat io at time t ^ 

Line 3: 

Field No. 

Field Name 

Format 

Columns 

Description 

1 

y 123 

E15.6 

2-16 

Ply 23 Poisson's ratio at time t ^ 

2 

l 'll3 

E15.6 

18-32 

Ply 13 Poisson's ratio at time t ^ 

3 

C fl 

E15.6 

34-48 

Ply heat capacity at time t ^ 

4 

K lll 

E15.6 

50-64 

Ply 11 thermal conductivity at 
time t. 

l 

5 

K 122 

E15.6 

66-80 

Ply 22 thermal conductivity at 
time t. 

i 
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Line 4: 


Field No. 

Field Name 

Format 

Columns 

1 

K 133 

E15.6 

2-16 

2 

“ill 

E15.6 

18-32 

3 

“l22 

E15.6 

34-48 

4 

“l33 

E15.6 

50-64 

5 

S 111T 

E15.6 

66-80 

Line 5: 




Field No. 

Field Name 

Format 

Columns 

1 

S 111C 

E15.6 

2-16 

2 

S 122T 

E15.6 

18-32 

3 

S 122C 

E15.6 

34-48 

4 

S 133T 

E15.6 

50-64 

5 

S 133C 

E15.6 

66-80 


Description 

Ply 33 thermal conductivity at 
time t^ 

Ply 11 coefficient of thermal 
expansion at time t^ 

Ply 22 coefficient of thermal 
expansion at time t^ 

Ply 33 coefficient of thermal 
expansion at time ^ 

Ply 11 tensile strength at time t^ 

Description 

Ply 11 compressive strength at 
t ime t 

Ply 22 tensile strength at time t^ 
Ply 22 compressive strength at 
time t. 

l 

Ply 33 tensile strength at time t. 
Ply 33 compressive strength at 
time t^ 
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Line 6: 


Field No. 

Field Name 

Format 

Columns 

Description 



1 

S 112 

E15.6 

2-16 

Ply 12 shear strength at 

time 

fc i 

2 

S 123 

E15.6 

18-32 

Ply 23 shear strength at 

time 


3 

S 113 

E15.6 

34-48 

Ply 13 shear strength at 

time 

l i 


Remarks : 

1. The set of six lines is repeated for each time step in the analysis. 
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4.6 Postprocessing File 6 


Description : Time and instantaneous fiber stresses and strains for each ply 


Format: 



Line 1: 


Field No. 

Field Name 

Format 

Columns 

Description 

i 


E15.6 

2-16 

Time at each load/time step 

2 

7 

E15.6 

18-32 

Number of fiber stresses 

3 

1 1 

E15.6 

34-48 

Fiber 11 stress at time t^ 

4 

°'f22 

E15.6 

50-64 

Fiber 22 stress at time t^ 

5 

*f\2 

E15.6 

66-80 

Fiber 12 shear stress at time t^ 
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Line 2: 


Field No. 

Field Name 

Format 

Columns 

Description 

1 

°f23 

E15.6 

2-16 

Fiber 23 shear stress at time t^ 

2 

°Yl3 

E15.6 

18-32 

Fiber 13 shear stress at time t^ 

3 

°f33 

E15.6 

34-48 

Fiber 33 stress at time t^ 

4 

7 

E15.6 

50-64 

Number of fiber strains 

5 

G fll 

E15.6 

66-80 

Fiber 11 strain at time t^ 


Line 3: 


Field No. 

Field Name 

Format 

Columns 

Description 

1 

e f 22 

E15.6 

2-16 

Fiber 22 strain at time t^ 

2 

G f 12 

E15.6 

18-32 

Fiber 12 shear strain at time t^ 

3 

C f23 

E15.6 

34-48 

Fiber 23 shear strain at time t^ 

4 

E f 13 

E15.6 

50-64 

Fiber 13 shear strain at time t^ 

5 

C f 33 

E15.6 

66-80 

Fiber 23 strain at time t^ 


Remarks : 

1. The set of three lines is repeated for each time step in the analysis. 
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4.7 Postprocessing File 7 


Description : Time and instantaneous matrix stresses and strains for each ply 


Format: 





+ 


+ 


+ 


+ 



19 


‘’mllA 


°'m22A 


°b22B 


* + 











"■22C 


°ml2A 


<r ml2B 


°'ml2C 




-T + 








_ — 


+ 

^n^B 


%23C 


^mlBA 


°ml3B 


^mlX 


T + 










+ 

^n^A 


%33B 


^330 


^mllB 


^rallC 


T i- + 










+ 

19 


C ml 1A 


G m22A 


C m22B 


C m22C 


+ 

C ml2A 


C ml2B 


e ml2C 


e m23A 


C m23B 

+ 

~ + 

C m23C 


C ml3A 


G ml3B 


G ml3C 


C m33A 

+ 

+ 
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G m33B 

C m33C 

e mllB 

-+ — — — — , 

e mllC 

Line 1 : 





Field No. 

Field Name 

Format 

Columns 

Description 

1 

l i 

E15.6 

2-16 

Time at each load/time step 

2 

19 

E15.6 

18-32 

Number of matrix stresses 

3 

^mllA 

E15.6 

34-48 

Matrix 11A stress at time tj 

4 

<7 'm22 A 

E15.6 

50-64 

Matrix 22A stress at time t 

5 

°m22B 

E15.6 

66-80 

Matrix 12B stress at time t^ 

Line 2: 





Field No. 

Field Name 

Format 

Columns 

Description 

1 

°m22C 

E15.6 

2-16 

Matrix 22B stress at time t^ 

2 

°ml2A 

E15.6 

18-32 

Matrix 12A shear stress at time t. 

3 

°ml2B 

E15.6 

34-48 

Matrix 12B shear stress at time t 

4 

% 12C 

E15.6 

50-64 

Matrix 12C shear stress at time t 

i 

5 

^2 3 A 

E15.6 

66-80 

Matrix 23A shear stress at time t 

i 

Line 3: 





Field No. 

Field Name 

Format 

Columns 

Description 

i 

°m23B 

E15.6 

2-16 

Matrix 23B stress at time t^ 

2 

^230 

E15.6 

18-32 

Matrix 23C shear stress at time t^ 

3 

‘’mOA 

E15.6 

34-48 

Matrix 13A shear stress at time 

4 

°ml3B 

E15.6 

50-64 

Matrix 13B shear stress at time t 

5 

^ml 3C 

E15.6 

66-80 

Matrix 13C shear stress at time 
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Line 4: 


Field No. 

Field Name 

Format 

Columns 

Descri 

ntion 




i 

°m33A 

E15.6 

2-16 

Matrix 

33A 

stress 

at 

time 


2 

°m33B 

E15.6 

18-32 

Matrix 

33B 

stress 

at 

time 

l i 

3 

^330 

E15.6 

34-48 

Matrix 

33C 

stress 

at 

time 

l i 

4 

^mllB 

E15.6 

50-64 

Matrix 

11B 

stress 

at 

time 


5 

^mllC 

E15.6 

66-80 

Matrix 

11C 

stress 

at 

time 

4 i 


Line 5: 


Field No. 

Field Name 

Format 

Columns 

Descri 

otion 




1 

19 

E15.6 

2-16 

Number 

of 1 

matrix j 

strains 


2 

e mllA 

E15.6 

18-32 

Matrix 

11A 

strain 

at 

time 

4 i 

3 

e m22A 

E15.6 

34-48 

Matrix 

22A 

strain 

at 

time 

fc i 

4 

e m22B 

E15.6 

50-64 

Matrix 

22B 

strain 

at 

time 

h 

5 

C m22C 

E15.6 

66-80 

Matrix 

22C 

strain 

at 

time 



Line 6: 


Field No. 

Field Name 

Format 

Columns 

Descri 

Dtion 





1 

C ml2A 

E15.6 

2-16 

Matrix 

12A 

shear 

strain 

at 

time 


2 

e ml2B 

E15.6 

18-32 

Matrix 

12B 

shear 

strain 

at 

time 

fc i 

3 

e ml2C 

E15.6 

34-48 

Matrix 

12C 

shear 

strain 

at 

time 

l i 

4 

C m23A 

E15.6 

50-64 

Matrix 

23A 

shear 

strain 

at 

time 

fc i 

5 

E m23B 

E15.6 

66-80 

Matrix 

23B 

shear 

strain 

at 

time 

fc i 
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Line 7: 


Field No 

Field Name 

Format 

Columns 

Description 

1 

G m23C 

E15.6 

2-16 

Matrix 23C shear strain at time t^ 

2 

e ml3A 

E15.6 

18-32 

Matrix 13A shear strain at time t^ 

3 

G ml3B 

E15.6 

34-48 

Matrix 13B shear strain at time t^ 

4 

e ml3C 

E15.6 

50-64 

Matrix 13C shear strain at time t^ 

5 

C m33A 

E15.6 

66-80 

Matrix 33A strain at time t^ 

Line 8: 





Field No. 

Field Name 

Format 

Columns 

Description 

1 

e m33B 

E15.6 

2-16 

Matrix 33B strain at time t^ 

2 

G m33C 

E15.6 

18-32 

Matrix 33C strain at time t^ 

3 

G mllB 

E15.6 

34-48 

Matrix 11B strain at time t^ 

4 

e mllC 

E15.6 

50-64 

Matrix 11C strain at time t. 

Remarks: 





1. The 

set of eight lines is 

repeated for 

each time step in the analysis. 
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4.8 Postprocessing File 8 


Description : Time and instantaneous interphase stresses and strains 

for each ply 


Format: 


n- i 

\ 

+ 

13 


°dllB 


°d22B 


^220 

' ■ < ~ 

<T 

dl2B 


^120 


°d23B 


°d23C 


°dl3B 

■ t -r “ — — 

( T 

d!3C 


°d33B 


°d33C 


°dllC 


13 


e dllB 


e d22B 


G d22C 


e dl2B 


e dl2C 

' T — 

e d23B 


e d23C 


e dl3B 


e dl3C 


e d33B 

. t » — — — — — — • — — -r — — — — — — — 

e d33C 


e dllC 








+ 


+ 


+ 


+ 


+ 
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Line 1 : 


Field No. 

Field Name 

Format 

Columns 

Description 


1 

fc i 

E15.6 

2-16 

Time at each load/time step 

2 

13 

E15.6 

18-32 

Number of interphase stresses 

3 

°dllB 

E15.6 

34-48 

Interphase 11B 

stress at time t^ 

4 

°d22B 

E15.6 

50-64 

Interphase 22B 

stress at time t^ 

5 

^220 

E15.6 

66-80 

Interphase 22C 

stress at time t^ 

Line 2: 

Field No. 

Field Name 

Format 

Columns 

Description 


1 

°dl2B 

E15.6 

2-16 

Interphase 12B 
time t^ 

shear stress at 

2 

°dl2C 

E15.6 

18-32 

Interphase 12C 
time t^ 

shear stress at 

3 

°d23B 

E15.6 

34-48 

Interphase 23B 
time t^ 

shear stress at 

4 

^230 

E15.6 

50-64 

Interphase 23C 
time t. 

l 

shear stress at 

5 

^dHB 

E15.6 

66-80 

Interphase 13B 
time t. 

l 

shear stress at 
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Line 3: 


Field No. 

Field Name 

Format 

Columns 

Descriotion 


1 

^130 

E15.6 

2-16 

Interphase 13C 

shear stress at 





time t^ 


2 

°d33B 

E15.6 

18-32 

Interphase 33B 

stress at time t 

3 

°d33C 

E15.6 

34-48 

Interphase 33C 

stress at time t 

4 

°dl 1C 

E15.6 

50-64 

Interphase 11C 

stress at time t 

5 

13 

E15.6 

66-80 

Number of interphase strains 

Line 4: 






Field No. 

Field Name 

Format 

Columns 

Description 


1 

e dllB 

E15.6 

2-16 

Interphase 11B 

strain at time t 

2 

e d22B 

E15.6 

18-32 

Interphase 22B 

strain at time t 

3 

C d22C 

E15.6 

34-48 

Interphase 22C 

strain at time t. 

4 

C dl2B 

E15.6 

50-64 

Interphase 12B 

shear strain at 





time t. 

l 


5 

e dl2C 

E15.6 

66-80 

Interphase 12C 

shear strain at 





time t^ 
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Line 5: 


Field No. 

Field Name 

Format 

Columns 

Descriotion 

1 

e d23B 

E15.6 

2-16 

Interphase 23B shear strain at 





time t^ 

2 

G d23C 

E15.6 

18-32 

Interphase 23C shear strain at 





time t^ 

3 

G dl3B 

E15.6 

34-48 

Interphase 13B shear strain at 





time t^ 

4 

C dl3C 

E15.6 

50-64 

Interphase 13C shear strain at 





time t. 

l 

5 

e d33B 

E15.6 

66-80 

Interphase 33B strain at time t^ 

Line 6: 





Field No. 

Field Name 

Format 

Columns 

DescriDtion 

1 

C d33C 

E15.6 

2-16 

Interphase 33C strain at time t 

2 

e dllC 

E15.6 

18-32 

Interphase 13C strain at time 

Remarks: 





1 . The set 

of six lines is repeated for 

each time step in the analysis. 
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4.9 Postprocessing File 9 


Description : Time, incremental composite strains, incremental ply stresses, 

incremental ply strains, and incremental composite stresses 


Format: 
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Line 1: 


Field No. 

Field Name 

Format 

Columns 

Description 


1 

t. 

i 

E15.6 

2-16 

Time at each load/time step 


2 

e 

cxx 

E15.6 

18-32 

Composite xx strain at time 


3 

e 

cyy 

E15.6 

34-48 

Composite yy strain at time 


4 

e 

cxy 

E15.6 

50-64 

Composite xy shear strain at 
time t^ 

5 

*111 

E15.6 

66-80 

Ply 11 stress at time t^ 


Line 2: 

Field No. 

Field Name 

Format 

Columns 

Description 


1 

*122 

E15.6 

2-16 

Ply 22 stress at time t^ 


2 

*112 

E15.6 

18-32 

Ply 12 shear stress at time 


3 

*113 

E15.6 

34-48 

Ply 13 shear stress at time 

t. 

l 

4 

*123 

E15.6 

50-64 

Ply 23 shear stress at time 


5 

c m 

E15.6 

66-80 

Ply 11 strain at time t^ 


Line 3: 

Field No. 

Field Name 

Format 

Columns 

Description 


1 

e l22 

E15.6 

2-16 

Ply 22 strain at time t^ 


2 

C 112 

E15.6 

18-32 

Ply 12 shear strain at time 

h 

3 

C 113 

E15.6 

34-48 

Ply 13 shear strain at time 

h 

4 

C 123 

E15.6 

50-64 

Ply 23 shear strain at time 

fc i 

5 

c r 

cxx 

E15.6 

66-80 

Composite xx stress at time 
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Line 4: 


Field No. 

Field Name Format 

Columns 

Description 

1 

O' 

cyy 

E15.6 

2-16 

Composite yy stress at time t 

2 

O' 

cxy 

E15.6 

18-32 

Composite xy shear stress at 





time t^ 

3 

O' 

CX2 

E15.6 

34-48 

Composite xz shear strain at 





time tj 

4 

O' 

cyz 

E15.6 

50-64 

Composite yz shear strain at 





time t^ 

5 

O' 

CZ2 

E15.6 

66-80 

Composite zz stress at time t^ 

Remarks: 





1. The set 

of four 

lines is repeated for 

each time step in the analysis. 
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4.10 Postprocessing File 10 


Description ; Time and instantaneous composite properties 


Format: 
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Line 1 : 


Field No. 

Field 

Name 

Format 

Columns 

Description 

1 

l i 


E15.6 

2-16 

Time at each load/time step 

2 

p c 


E15.6 

18-32 

Composite density at time t^ 

3 

T 

c 


E15.6 

34-48 

Composite temperature at time 

4 

C cll 


E15.6 

50-64 

3-D composite 11 component of 
stiffness matrix at time t^ 

5 

C cl2 


E15.6 

66-80 

3-D composite 12 component of 
stiffness matrix at time t^ 

Line 2: 

Field No. 

Field 

Name 

Format 

Columns 

Description 

i 

C cl3 


E15.6 

2-16 

3-D composite 13 component of 
stiffness matrix at time t^ 

2 

C c22 


E15.6 

18-32 

3-D composite 22 component of 
stiffness matrix at time t^ 

3 

C c23 


E15.6 

34-48 

3-D composite 23 component of 
stiffness matrix at time t^ 

4 

C c33 


E15.6 

50-64 

3-D composite 33 component of 
stiffness matrix at time t. 

l 

5 

C c44 


E15.6 

66-80 

3-D composite 44 component of 
stiffness matrix at time t^ 
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Line 3: 


Field No. 

Field Name 

Format 

Columns 

1 

C c55 

E15.6 

2-16 

2 

C c66 

E15.6 

18-32 

3 

“ell 

E15.6 

34-48 

4 

a c22 

E15.6 

50-64 

5 

a c33 

E15.6 

66-80 

Line 4: 




Field No. 

Field Name 

Format 

Columns 

1 

K cll 

E15.6 

2-16 

2 

K c22 

E15.6 

18-32 

3 

K c33 

E15.6 

34-48 

4 

C fc 

E15.6 

50-64 

5 

E cll 

E15.6 

66-80 


Description 

3-D composite 55 component of 

stiffness matrix at time t. 

1 

3-D composite 66 component of 
stiffness matrix at time t^ 

3-D composite 11 coefficient of 
thermal expansion at time t^ 

3-D composite 22 coefficient of 
thermal expansion at time t^ 

3-D composite 33 coefficient of 
thermal expansion at time t^ 

Description 

3-D composite 11 thermal 
conductivity at time t^ 

3-D composite 22 thermal 
conductivity at time t^ 

3-D composite 33 thermal 
conductivity at time t^ 

3-D composite heat capacity at 
time t^ 

3-D composite 11 modulus at time t^ 
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Line 5: 


Field No. 

Field Name 

Format 

Columns 

1 

E c22 

E15.6 

2-16 

2 

E c33 

E15.6 

18-32 

3 

G c23 

E15.6 

34-48 

4 

G c31 

E15.6 

50-64 

5 

G c12 

E15.6 

66-80 

Line 6: 




Field No. 

Field Name 

Format 

Columns 

1 

U cl2 

E15.6 

2-16 

2 

V c 2 l 

E15.6 

18-32 

3 

U cl3 

E15.6 

34-48 

4 

V c31 

E15.6 

50-64 

5 

V c23 

E15.6 

66-80 


Description 

3-D composite 22 modulus at time t. 

1 

3-D composite 33 modulus at time t 

l 

3-D composite 23 shear modulus at 
time t^ 

3-D composite 31 shear modulus at 
time t^ 

3-D composite 12 shear modulus at 
time t. 

l 

Description 

3-D composite 12 Poisson’s ratio 
at time t^ 

3-D composite 21 Poisson’s ratio 
at time t^ 

3-D composite 13 Poisson’s ratio 
at time t. 

l 

3-D composite 31 Poisson’s ratio 
at time t. 

l 

3-D composite 23 Poisson’s ratio 
at time t^ 
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Line 7: 


Field No . Field Name Format Columns Description 

1 y „ E15.6 2-16 3-D composite 32 Poisson’s ratio 

c32 

at time t. 

l 

2 z c E15.6 18-32 Composite thickness 

3 \b E15.6 34-48 Composite modified distortion 

*cMDE 

energy 

4 C c n E15.6 50-64 2-D composite 11 component of 

stiffness matrix at time t. 

l 

5 ^cl2 E15.6 66-80 2-D composite 12 component of 

stiffness matrix at time t^ 

Line 8: 

Field No. Field Name Format Columns Description 

1 ^cl3 E15.6 2-16 2-D composite 13 component of 

stiffness matrix at time 

2 C c22 E15.6 18-32 2-D composite 22 component of 

stiffness matrix at time ti 

3 ^c23 E15.6 34-48 2-D composite 23 component of 

stiffness matrix at time t. 

i 

4 ^c33 E15.6 50-64 2-D composite 33 component of 

stiffness matrix at time t^ 

5 E E15.6 66-80 2-D composite 11 modulus at time i. 

cl 1 i 
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Line 9: 


Field No. 

Field Name 

Format 

Columns 

1 

E c22 

E15.6 

2-16 

2 

G c12 

E15.6 

18-32 

3 

V cl2 

E15.6 

34-48 

4 

V c21 

E15.6 

50-64 

5 

l 'cl3 

E15.6 

66-80 

-ine 10: 

■ield.No. 

Field Name 

Format 

Columns 

1 

l c31 

E15.6 

2-16 

2 

V c23 

E15.6 

18-32 

3 

U c32 

E15.6 

34-48 

4 

“ell 

E15.6 

50-64 

5 

a c22 

E15.6 

66-80 


Description 

2-D composite 22 modulus at time t^ 
2-D composite 12 shear modulus at 
time t^ 

2-D composite 12 Poisson's ratio 
at time t. 

l 

2-D composite 21 Poisson’s ratio 
at time t^ 

2-D composite 13 Poisson’s ratio 
at time t. 

l 

Description 

2-D composite 31 Poisson’s ratio 
at time t^ 

2-D composite 23 Poisson’s ratio 
at time t. 

l 

2-D composite 32 Poisson’s ratio 
at time t. 

l 

2-D composite 11 coefficient of 
thermal expansion at time t 

i 

2-D composite 22 coefficient of 
thermal expansion at time t 

i 
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Line 11: 


Field No. 

Field Name 

Format 

Columns 

Description 

1 

°cl2 

E15.6 

2-16 

2-D composite 12 coefficient 
thermal expansion at time t^ 

2 

K cll 

E15.6 

18-32 

2-D composite 11 thermal 
conductivity at time t^ 

3 

K c22 

E15.6 

34-48 

2-D composite 22 thermal 
conductivity at time t^ 

4 

K c12 

E15.6 

50-64 

2-D composite 12 thermal 
conductivity at time t^ 

5 

C fc 

E15.6 

66-80 

2-D composite heat capacity 


Line 12: 


Field No. 

Field Name 

Format 

Columns 

Description 



1 

vacant 

E15.6 

2-16 

reserved 

for 

future 

use 

2 

vacant 

E15.6 

18-32 

reserved 

for 

future 

use 

3 

vacant 

E15.6 

34-48 

reserved 

for 

future 

use 

4 

vacant 

E15.6 

50-64 

reserved 

for 

future 

use 

5 

vacant 

E15.6 

66-80 

reserved 

for 

future 

use 


Line 13: 


Field No. 

Field Name 

Format 

Columns 

Description 



1 

vacant 

E15.6 

2-16 

reserved 

for 

future 

use 

2 

vacant 

E15.6 

18-32 

reserved 

for 

future 

use 

3 

vacant 

E15.6 

34-48 

reserved 

for 

future 

use 


Remarks : 

1. The set of thirteen lines is repeated for each time step in the analysis. 
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